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Abstract 



We analyse three different New Physics scenarios for AF = 2 flavour-changing neutral 
currents in the quark sector in the light of recent data on neutral-meson mixing. We 
parametrise generic New Physics contributions to B q — B q mixing, q — d,s, in terms of 
one complex quantity A q , while three parameters A^-, and A c £ are needed to describe 
K — K mixing. In Scenario I, we consider uncorrelated New Physics contributions in the 
B d , B s , and K sectors. In this scenario, it is only possible to constrain the parameters A d 
and A s whereas there are no non-trivial constraints on the kaon parameters. In Scenario 
II, we study the case of Minimal Flavour Violation (MFV) and small bottom Yukawa 
coupling, where A = A^ = A s = A^. We show that A must then be real, so that no new 
CP phases can be accomodated, and express the remaining parameters A^ and A^ in 
terms of A in this scenario. Scenario III is the generic MFV case with large bottom Yukawa 
couplings. In this case, the Kaon sector is uncorrelated to the Bd and B s sectors. As in 
the second scenario one has Ad = A s = A, however, now with a complex parameter A. 
Our quantitative analyses consist of global CKM fits within the Rfit frequentist statistical 
approach, determining the Standard Model parameters and the new physics parameters 
of the studied scenarios simultaneously. We find that the recent measurements indicating 
discrepancies with the Standard Model are well accomodated in Scenarios I and III with 
new mixing phases, with a slight preference for Scenario I that permits different new CP 
phases in the Bd and B s systems. Within our statistical framework, we find evidence of 
New Physics in both B d and B s systems. The Standard-Model hypothesis A d = A s = 1 
is disfavoured with p-values of 3.6(7 and 3.3a in Scenarios I and III, respectively. We 
also present an exhaustive list of numerical predictions in each scenario. In particular, we 
predict the CP phase in B s — > J/ip4> and the difference between the B s and Bd semileptonic 
asymmetries, which will be both measured by the LHCb experiment. 
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1 Introduction 



1 Introduction 

Considerations of the stability of the electroweak scale lead to the general belief that 
there is New Physics with particle masses below 1 TeV. While the high-p^ experiments at 
the LHC should produce these new particles directly, one can study their dynamics also 
indirectly, through their impact on precision measurements at lower energies. To this end 
flavour-changing neutral current (FCNC) processes are extremely useful. On one hand 
they are highly suppressed in the Standard Model and are therefore very sensitive to new 
physics. On the other hand FCNC processes of K, and B s mesons are still large enough 
to be studied with high statistics in dedicated experiments. Here meson-ant imeson mixing 
plays an outstanding role. First, meson-antimeson oscillations occur at time scales which 
are sufficiently close to the meson lifetimes to permit their experimental investigation. 
Second, the Standard Model contribution to meson-antimeson mixing is loop-suppressed 
and comes with two or more small elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix pQ. Third, the decays of oscillating mesons give access to many mixing- induced 
CP asymmetries through the time- dependent study of decays into CP-eigenstates, which 
in some cases one can relate to the parameters of the underlying theory with negligible 
hadronic uncertainties. 

The 5-factories have revealed that the dominant b — > d and b — > u transitions fit 
into the pattern of the CKM mechanism and are in agreement with the information on 
s — > d transitions gained in more than fourty years of kaon physics. The success of the 
CKM picture is evident from the many different measurements combining into a consistent 
and precise determination of the apex (p, rj) of the S-meson unitarity triangle (in terms 
of the Wolfenstein parameterization of the CKM matrix [21 E]). As a consequence, any 
contribution from the expected new TeV-scale physics to the measured flavour- changing 
processes must be suppressed compared to the established CKM mechanism. 

Models with only CKM-like flavour violation are said to respect the principle of min- 
imal flavour violation (MFV) [UG3. This principle is often invoked in an ad- hoc way to 
suppress excessive FCNC amplitudes for model-building purposes. In the Minimal Su- 
persymmetric Standard Model (MSSM), new sources of flavour violation solely stem from 
the supersymmetry-breaking sector. A sufficient condition for MFV are supersymmetry- 
breaking terms which are flavour-blind at a given energy scale. This situation occurs in 
supergravity with a flat Kahler metric [B] or if supersymmetry breaking is mediated by 
gauge interactions [7] . The overall picture of experimental data does not require sizeable 
corrections to MFV. Still it is difficult to probe the CKM picture with a better accuracy 
than, say, 30%, because most quantities entering the global fit of the unitarity triangle 
suffer from sizeable hadronic uncertainties. It should also be stressed that the accuracy of 
the determination of the CKM parameters decreases notably when one assumes that one 
or several crucial input (s) could be affected by New Physics contributions. Interestingly, 
several authors have detected possible hints of New Physics in the data. For example it 
has been argued in the literature that one starts to see a discrepancy between the mea- 
surement of sin 2/3 and the region preferred by \V u b\ from semileptonic decays on one hand, 
and \£k\ on the other hand [8,9J. Also the recently improved measurement of the B — )■ rv 
branching ratio deviates from its indirect CKM fit prediction [TO]. In addition there are 
anomalies in the data on b — > s transitions. The latter processes do not involve p and fj (to 
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a good accuracy) and therefore directly probe the CKM mechanism. An ongoing debate 
addresses an extra contribution to b — > sqq, q = u,d, s, decay amplitudes with a CP phase 
different from aig(V t * s V tb ) that can alleviate the pattern of shifts between the measured 
CP asymmetries in these b — > s penguin modes and the Standard Model predictions (see 
e.g. (HI). 

However, the first place to look for New Physics in b — > s transitions is B s —B s mixing, 
where New Physics can be parameterised in terms of just one complex parameter (or two 
real parameters) in a mo del- independent way, as we will discuss in great detail below. At 
the end of 2006 a combined analysis of several observables has pointed to the possibility of 
a new-physics contribution with a CP phase different from that of the Standard-Model box 
diagram [12]. Models of supersymmetric grand unification can naturally accommodate new 
contributions to b — > s transitions [13] : right-handed quarks reside in the same quintuplets 
of SU(5) as left-handed neutrinos, so that the large atmospheric neutrino mixing angle 
could well affect squark-gluino mediated b — > s transitions. At the same time the GUT 
models of Refs. [T3] do not induce too dangerous contributions to the well-measured rare 
decay B — > X s ^. B s — B s mixing has been further investigated in other supersymmetric 
scenarios with [TJ] and without [15j[T6] GUT boundary conditions, in unparticle physics 
scenarios in Ref. [17], in multi-Higgs-doublet models [18], in models with extra gauge 
bosons Z' [19], warped extra dimensions [20] , left-right symmetry [21], anomalous tWb- 
couplings [22], additional quark families [23] or an additional singlet quark [21], and in a 
little-Higgs model [26]. 

On the experimental side, the understanding of b — > s transitions has made a tremen- 
dous progress in the past years. The Tevatron experiments have discovered and precisely 
quantified B s — B s mixing oscillations [27|l2"8] whose frequency is in good agreement with 
the Standard Model prediction, and presented first determinations of the associated CP- 
violating phase from tagged analyses of B s — > J/ip4> decays [2"9"H3"3] . Recently, possible 
New Physics in the CP phase of the B s — B s mixing amplitude has received new atten- 
tion: The D0 collaboration has reported a measurement of the dimuon charge asymmetry 
which disagrees with the Standard Model prediction by 3.2 standard deviations |35j. The 
CP asymmetry in semileptonic or, more generally, any flavour-specific decays, involves the 
B s —B s mixing phase just as B s — > J/ip4>, so that both pieces of experimental information 
can be combined to constrain this phase. The new measurement of the dimuon charge 
asymmetry has already triggered considerable theoretical interest. Besides predictions for 
the CP phase of B s —B s mixing in specific models, as quoted above, also model-independent 
analyses of New Physics effects have appeared [36]. Due to the large size of the dimuon 
asymmetry it was also investigated whether sizeable New Physics contributions to the 
decay of B s mesons are possible. This alternative is however strongly constrained by the 
lifetime ratios of B- mesons (see e.g. [37]) as well as the semileptonic branching ratios and 
the the average number of charm quarks per b-decays (see e.g. [38]). 

In the present article, we analyse generic scenarios of New Physics which are compatible 
(at different levels) with the experimental picture sketched above. We set up our notation 
and define our theoretical framework in Sect. [2j the relevant updated experimental and 
theoretical inputs to our global analysis are presented in Sect. |3J In Sect. 14.11 we first 
present the current status of the Standard Model fit. In Sect. 14.21 we perform a fit in 
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which we allow for New Physics in the Bd — Bd mixing and the B s — B s mixing systems 
and we project the results onto the New Physics parameters that describe the Bd — Bd 
mixing and the B s — B s mixing systems. Sects. 14.31 and 14.41 are dedicated to two MFV 
scenarios with correlated effects in all meson-antimeson mixing amplitudes. Finally, we 
conclude and list a few perspectives for the close future. 

2 Setting the scene 
2.1 B — B mixing basics 

B q — B q oscillations (with q = d or q = s) are described by a Schrodinger equation: 



:! ( \B&)) \_ ( Mq _ l TQ \ ( \B&)) 

dt I \B q (t)) ) ~ I 2 1 J I \B q (t)) 



(1) 



with the mass matrix M q = M q ^ and the decay matrix T q = T q \ The physical eigenstates 
\B q H ) and \Bf) with masses Mjj, Ml and decay rates T q H , T q L are obtained by diagonalizing 
M q — iT q /2. The B q —B q oscillations in Eq. ([T|) involve the three physical quantities |Mf 2 |, 
\V\ 2 \ and the CP phase 

9 = arg(-Mf 2 /r? 2 ). (2) 

We denote the average B q mass and width by M Bq and T Bq , respectively. The mass and 
width differences between B q L and B q H are related to them as 

AM q = M q H -M q L = 2\M q 2 \, ATq = T q L — T q H = 2 |rf 2 | COS0,, (3) 

up to numerically irrelevant corrections of order m^/M^. AM, simply equals the fre- 
quency of the Bq — Bq oscillations (for details see e.g. [39]). A third quantity probing 
mixing is 



4 = 2(1 



T q \T q I Ar 
I m vr^r — i , la i sin d) g = - q tan0 o . (4) 

M 12 \ M 12\ AM q 



a q s is the CP asymmetry in flavour-specific B q — )■ / decays, i.e., the decays B q —y f and 
B q f are forbidden. The standard way to measure a q s uses B q — > Xt+v^ decays, which 
explains the common name semileptonic CP asymmetry for a q s , with the corresponding 
notation ag L (for more details see e.g. [40J). In theoretical contexts, we use the notation 
aj s in this paper, while we write a| L when referring to the specific experimental observable 
inferred from semileptonic decays. Further 



q 2M? 2 -iT* 12 



p V 2 M12 - iT 



(5) 



12 



Let us now discuss our theoretical understanding of the off-diagonal terms of the evolution 
matrix, which are responsible for B q — B q mixing. The dispersive term M\ 2 is completely 
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dominated by box diagrams involving virtual top quarks, and it is related to the effective 
| AS | = 2 Hamiltonian ifj, AB ' _2 as 

The Standard Model expression for H l q AB] -' IS 1411 

hI abi = 2 = (y; q v tb fcQ + h.c. (7) 

with the four-quark operator 

Q = q L l^Lq L Y b L, (8) 
q L = \{l-l,)q, (9) 

and the Wilson coefficient C, which depends on the heavy mass scales of the theory. In 
a wider class of models .£/j AB ' -2 maintains the form of Eq. (j7j) (meaning that there is no 
other operator than Q involved), but with a value of C different from the value in the 
Standard Model: 

Gl . , _ „ / m. 2 
Mr 



C- = - ). (10) 



Here m f is the top quark mass defined in the MS scheme, related to the pole mass mj? oIe 
determined at the Tevatron as m t (jn t ) = 0.957 m]? ole (at next-to-leading order of QCD). 
The Inami-Lim function S [42J is calculated from the box diagram with two internal top 
quarks and evaluates to Simf/M^) = 2.35 for the central value of m t listed in Table [71 
QCD corrections are comprised in [I3j| 



rj B = 0.8393 ± 0.0034. (11) 
The hadronic matrix element involved is usually parameterised as 

(B q \Q^B)\B q ) = \Ml q f B B Bq {VB), (12) 

with the decay constant f Bq and the 'bag' factor B Bq - The product rj B B Bq is scale and 
scheme invariant. Our convention in Eq. ffTTl) corresponds to a scale dependent bag pa- 
rameter with B B = 1 in vacuum insertion approximation. Typical values for the bag 
parameter obtained on the lattice are e.g. B Bs ~ 0.84, see Sect. 13.11 Sometimes a different 
normalisation with a scale independent bag parameter B Bq = b B (fi B )B Bq (fi B ) is used. The 
corresponding quantities r\ B = rj B (^ B )/b B (fj, B ) = 0.551 and B Bs ~ 1.28 obviously satisfy 
VbB Bs = rj B B Bs . The analytic formula for b B (fi B ) can be found e.g. in Eq. (XIII. 5) of [45] . 

The absorptive term T\ 2 is dominated by on-shell charmed intermediate states, and 
it can be expressed as a two-point correlator of the |A5| = 1 Hamiltonian i^ 5 ' -1 . By 
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performing a 1/ m^-expansion of this two-point correlator, one can express rf 2 in terms of 
Q and another four-quark operator 

Qs = q a L b R f L b a R , (13) 

where S stands for "scalar" and a, (3 = 1,2,3 are colour indices, see [12]. The matrix 
element is expressed as 

(B,\Qs\B,) = ^^(^y^ij^/^. (14) 

The prediction of T\ 2 involves also operators which are subleading in the heavy quark 
expansion, the matrix elements of which are parameterised by the 'bag' factors B Roi23 
and B Ri [12] , which satisfy two relations in the heavy quark limit [12j[68] : 

B R2 =B R2 , B R3 = ^-B R3 + 2 -B k2 . (15) 

Even though we have not included the flavour of the light-quark in our notation, we 
consider these l/m^-suppressed operators to have different values for B d and B s mesons. 

Finally, we discuss the relative phase <p q between the two off-diagonal terms. In contrast 
to M? 2 , Y\ 2 receives non-negligible contributions from subleading u and c CKM couplings, 
which implies that 4> q is not a pure CKM phase in the Standard Model. The Standard 
Model contribution to 4> q reads [121116], with our updated inputs (see Table [9]) 

= (-io.iiL 7 )-io- 2 , 

SM = (+7.4+0.8) . 1Q -3 ; (16) 

and thus in <fr s the Standard Model contribution is clearly subleading in the presence of 
generic New Physics effects. 

The previous quantities are expected to be affected by New Physics in different ways. 
While Mf 2 coming from box diagrams is very sensitive to New Physics both for Bd and 
B s , T{ 2 stems from Cabibbo-favoured tree-level decays and possible New Physics effects 
are expected to be smaller than the hadronic uncertainties. In the case of Tf 2 though, 
the contributing decays are Cabibbo-suppressed. In this paper, we only consider scenarios 
where New Physics does not enter tree-level decays. More specifically, we assume that B 
decays proceeding through a four-flavour change (i.e., b — >■ qi&qs, Qi 7^ <l2 ^ <h) obtain 
only Standard Model contributions (SM^FC) [171HE] . This assumption is better defined 
than just the neglect of New Physics contributions to tree-mediated decays since on the 
non-perturbative level tree and penguin amplitudes can not be well separated. Our class 
of four-flavour-change decays includes b — > d decays in which the strong isospin changes 
by 3/2 units, i.e. we use strong isospin as the flavour quantum number of the first quark 
generation. Then the following inputs used in the fit are considered to be free from New 
Physics contributions in their extraction from data: \V u d\, \V US \, \V U b\, \V c b\ and 7. Also 
the leptonic decays B — > tu (or D s — > rv and D s — > /iz/), which could be significantly 
affected by charged Higgs exchange contributions (see [1H] and references within), are 
assumed to be Standard Model-like. Using these inputs a reference unitarity triangle 
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can be constructed (see the first two articles in Ref. [17]). as will be discussed further in 
Sect. 13.21 (in Ref. [50], this triangle is compared with the universal unitarity triangle for 
models of constrained MFV introduced in Ref. [5T]). 

In addition, in order to take advantage of the measurement of the width differences 
Ar g , and of the time-dependent CP-asymmetries in dominant b — > c decays, we neglect 
possible non standard contributions to the b — > ccq (q = d, s) transitions, although they 
do not strictly enter the SM4FC family. Finally, we assume that the unitarity of the 
3 x 3-CKM matrix still holds in the presence of New Physics, which ensures that the 
Standard Model contribution to the neutral meson mixing keeps its usual expression as a 
function of (p, fj) and other parameters. Hence, our discussion would not hold in the case 
of an additional sequential fourth fermion family, which however is not excluded yet by 
experimental constraints (see Refs. [23J or [52J and references therein). 

Thus, New Physics can find its way into the quantities studied in this paper only by 
changing magnitude and/or phase of Mf 2 . It is convenient to define the New Physics 
complex parameters A q and (fi^ (q = d, s) through 



M\ 2 ee M™' q -A q , A q = \A q \e^, (17) 

see e.g. [32] • With the definition in Eq. ( [171) the CP phase of Eq. ([2|) reads 

0, = € M + K- ( 18 ) 

As discussed in Sect. [31 New Physics in M\ 2 will not only affect the neutral- meson mixing 
parameters, but also the time-dependent analyses of decays corresponding to an interfer- 
ence between mixing and decay. 

The relation of A q to the parameters used e.g. in [531455] is | A 9 | = r q , cj) q = 29 q , and 
the Standard Model is of course located at A q = 1. It is more transparent to look at 
the cartesian ImA, vs. Re A^ plot than the polar 26 q vs. r q one, because it visualizes the 
New Physics contribution more clearly and it allows a simple geometrical interpretation 
of the shape of each individual constraint. For completeness, we note that some authors 
(e.g. [2BHSE], see a l so [S]) prefer to split the Standard Model contribution from the pure 
New Physics one in a polar parametrization. The two New Physics parameters h q and 2a q 
introduced in this way are defined by 

ML M™'" A :n_ 

-i + dtr^ A, = l + h q e^ . (19) 



Mf 2 M ' 9 M*F* q ~ q 

We will study the case of the neutral kaon system in the following section, defining ana- 
logous parameters A^-, A^ and A^. But in this paper, we will not consider the neutral 
D meson system. Indeed, in the scenarios we consider here D — D mixing is severely 
GIM-suppressed and gives no useful constraint, as is already the case within the Standard 
Model, see e.g. 



2.2 K — K mixing basics 

The effective \AS\ = 2 Hamiltonian describing K — K mixing resembles the \AB\ = 2 
Hamiltonian of Eq. (jTJ, with the important distinction that now also contributions from 
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internal charm quarks are important: 

h \as\=2 = [(VuV^Cu + 2v ts v t * d v cs v; d c ct + (v cs v; d ) 2 c cc ] q + h.c. (20) 

with the operator Q = cLl^^sl cLl^sl 0- As for the case of the B d and B s mesons, the 
contribution from ifl A5 l= 2 to Mf^ is found from Eq. ([6]). A new feature is an additive 
poorly calculable long-distance contribution involving ifl AS 1 =1 (see e.g. ED] ) - The 
Wilson coefficients CV,, i,j = c, t, and the operator Q depend on the renormalisation scale 
fiK at which we evaluate these coefficients and the hadronic matrix element (K\Q\K). We 
parametrise the hadronic matrix element as 

(K\QM\K) = \ m \f K -^- (21) 

Here fx — 156 MeV and mx are the decay constant and mass of the kaon, respectively, 
and Bk is the bag parameter, from which a factor ^(//r-) is stripped off. Analogous to the 
case of B mixing bxi^x) contains the dependence of (K\Q(/j,k)\K) on the renormalisation 
scheme and the renormalisation scale fix- The Standard Model values of the Wilson 
coefficients are 

= ^M> w s(^^-) Vct b K (» K ), 



4vr 2 w \M^M W 
<s\j G 2 F 2 „ ( m 2 



C b cc M = ^M 2 ,S[ J ^) Vcc b K ( f i K ), (22) 



with the Inami-Lim functions S calculated from the usual box diagrams. By comparing 
Eqs. ([IUD and §M) we verify the MFV feature of the Standard Model C SM = Cf t M . In C C S C M 
the Inami-Lim function can be expanded in terms of the tiny quantity m 2 /My / to find 
S(x c ) = x c + 0(x 2 ). Likewise S(x c ,x t ) ~ — x c logx c + x c F(x t ) with F(m 2 /M^) = 0.56. 
From Eqs. f l2T]) and ff22|) . bjcifiR:) drops out if (K\H^ AS ^ =2 \K) is expressed in terms of 
Bk- The QCD correction factors r\ cc [61], r\ ct [62] and r] tt [13] are listed in Table [3 The 
dominant sources of uncertainty in these quantities are higher-order QCD corrections^] (i] cc 
also depends on a s and fn c in a sizeable way). The latter dependence is made explicit in 
Tabled 

In analogy to Eq. (TPTj) we introduce complex parameters for New Physics in the three 
different contributions and write 

m« = [ 1 2Mk 1 ; = {v ts v: d ) 2 M^+2v ts v: d v cs v: d Mti + {v cs v: d fMi^ 

MH = Afg^'Ag = Mff^lAlle 1 ^. (23) 



: We use the same notation for operators in the B and K systems (cf. Eq. (|8j)). implying the corre- 
sponding flavours (b, s or d) of the quark fields. 

2 Very recently a NNLO calculation of rj c t was performed [63], leading to a value, which is 5% larger 
than the value used here. This result is not yet included in our analysis. 
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The physical quantities associated with K — K mixing are the K^-Ks mass difference 
AM K = M Kl - M Ks and the CP-violating quantity e K . CP violation in ifl A5 l =1 is 
characterised by e' K . These quantities are defined as 



ex 



r] 00 + 2rj + _ 



-Voo + V+- 



A[K L -> tt VQ 
A(K S ->■ n a ii b ) ' 



with T), 



ab 



(24) 



Since these two quantities are defined in terms of Kl and Ks, they can be expressed in 
terms of K — K mixing parameters and the isospin decay amplitudes A(Kq — > = 
Aje tSl = aie t6i e lSl , where a/, 5j and 9j denote the modulus, the "strong" (CP-even) phase 
and the "weak" (CP-odd) phase of the decay amplitude [391 EH ES] . e# is essentially 
proportional to the CP phase 4>k = ar g( — ^12/^12) ■ ^ n y i ew °f the phenomenological 
"A/ = 1/2 rule" a /a 2 ~ 22 (and the fact that all other decay modes come with even 
smaller amplitudes than 02) one can saturate the inclusive quantity T{| completely by the 
contribution proportional to a^. Expanding in various small parameters (see Ref. [39] for 
an elaborate discussion of the approximations involved) one finds: 



sin 1 



AM, 



K 



with tan < 



2AM K 
AT 



and £ 



Im A, 







A 



Re A r 



(25) 



The troublesome long-distance contribution to mentioned after Eq. ( 120]) is eliminated 
from Eq. (|25|) by trading 2ReM^ for the experimental value of AM^. Long-distance 
contributions to ImMj| are small |66j. In Eq. (125]) £ comprises the contribution from 
arg(— T^) in the limit of A dominance discussed above. The corrections are of order 
(o^/ao) 2 an d therefore negligible. The usual expression for ex is obtained from this ex- 
pression by taking the following further approximations: i) use <p e = 45° instead of the 
measured value 4> e = 43.5(7)°, ii) neglect £ and iii) compute ImM 12 using only the lowest- 
dimension d = 6 operator in the effective Hamiltonian of Eq. ( 120]) . which is dominated by 
top and charm box diagrams. The effect of the three simplifications can be parameterised 
in terms of the parameter n e [9] entering 



ImM 



(6) 

12 



AM 



(26) 



C e K e e^B K 



Im [(V cs V; d ) 2 A-] Vcc S 



— 2 

m c 



+ Im [{V ts V; d Y A«] Vtt S 



— 2 

m t 



+ 2Im (V ts V t * d V cs V* d A c x) Vct S 



2 2 



lvl w w 



(27) 



The value K e = 1 corresponds to the approximations i) — iii) outlined above. The normal- 
isation reads 



a 



G 2 F F 2 K m K M^ 



{2i 



12V2tt 2 AM k ' 

When expressed in terms of Wolfenstein parameters to lowest order in A, Eq. ([26]) defines 
the familiar hyperbola in the ~p-rj plane. 

A series of papers [91I31I ] I6"6" ] I6"7] has studied how much the factor n e should deviate from 
1 in order to account for the terms neglected by the previous approximations. We recall 
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the different elements in App. |XJ separating the uncertainties coming from statistical and 
systematic sources, and we obtain the estimate 

n t = 0.940 ±0.013 ±0.023, (29) 

in good agreement with K t = 0.94 ± 0.02 in Ref. [66]. We emphasize that the estimate of 
k 6 in Eq. ( 1291) relies on the assumption that e' K is unaffected by New Physics (which goes 
beyond the SM4FC assumption which protects only 1 = 2 final states). From ex,exp = 
(2.229 ± 0.010) ■ 10~ 3 we get the following value for the combination: 

eP exp = = (2.367 ± 0.033 ± 0.049) • KT 3 . (30) 

In the presence of New Physics, the relationship between the measured ex and the A^-'s 
is discussed after Eq. (1391 . 

One can also study the semileptonic CP asymmetry 

. _ T(K long ^e+vn-)-T(K long ^e-V7r + ) 1 - \q/p\ 2 

At 



T(K long -> + T(K long -> £-un+) 1 + |g/p| 2 ' 

which, however, contains the same information on fundamental parameters as Ree^. 

2.3 Master formulae 

In this section, we provide the master formulae of the theoretical predictions for the 
observables relevant to the analysis of new-physics contributions in mixing. These formulae 
reflect the dependences on the most important parameters entering the fit and are obtained 
from the input values as described in Sect. [3j It should be stressed that these numerical 
equations are shown for illustrative purpose only: the complete formulae are used in the 
fitting code, which allow to take into account all the contributions computed so far together 
with the correct treatment of the correlations. 

Combining Eqs. fl3T 4T2|) with Eq. ( ITT)) one finds: 

1 V 0.0086 7 2.35 (0.17 GeV) 2 1 

AM = 17 24 P s- 1 YI^V gM^ ^ -|A| (31) 
' PS V 0.04 J 2.35 (0.21 GeV) 2 1 s| ' 1 j 

The remaining uncertainties in the prefactors of the above formulae are due to the choice 
of the renormalization scale and the values of a s and the top quark mass. They are at most 
3 % and therefore negligible compared to the theoretical error due to the non-perturbative 
and CKM parameters. 

The derivation of the formulae involving T\ 2 is more complicated [T2Tl46y 68] . For the B s 
system the dependence on the apex (p, fj) of the unitarity triangle is strongly suppressed, 
in contrast to the B d system. Furthermore the Standard Model contribution to af s is tiny 
and remains below the present experimental sensitivity, while af s is one order of magnitude 
larger and therefore not completely negligible, see Table O Summing up logarithms of the 
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form m 2 c /m1 him 2 /ml [69] and using the MS-scheme for the b quark mass one finds from 
Ref. [12] for the decay rate differences: 



at. 



I 0.17 GeV / 



B' B 
0.00241 + 0.00056^+^ - 0.00047- R,Bd 



B 



B d 



cos 



M M + d A ) , (32) 



0.21 GeV 



B' 



B 

0.0797 + 0.0278^^ - 0.0181- VM ' 



B 



B 



B s 



cos (<^ M + <£) . (33) 



Now the uncertainties in the coefficients are considerably larger than in the case of the 
mass differences, but they are still less than about 15 %. The dominant theoretical error of 
the coefficients comes from the renormalization scale pi followed by the CKM factors. One 
encounters matrix elements of higher-dimensional operators in these expressions, denoted 
by Br, which have a power suppression parametrised by m v ™ . The general (assuming 
unitarity of the 3x3 quark mixing matrix) expression for the semileptonic CP asymmetries 
reads 



io 4 < 



a g Im 



+ b a \m 



a; 7 



sin 



with X q x = V x i)V* . The coefficients a, b, c read 



(34) 



b d 
b s 

Cd 



9.2905 + 0.2973- 
9.4432 + 0.2904 
0.0720 + 0.0184 
0.0732 + 0.0180- 



5 ^ + 0.2830 



Be 



Bra 
B Bd 



B{ 



S,B S , r. nrrnBR,B s 



Be 



B 



+ 0.2650 
+ 0.0408 



B Bs 
Bn,B d 



s,b ± 

B d VB d 



S,B S 



Bi 



+ 0.0395 



Bi 

Br,b s 



-46.8169 - 17.0083- 



B, 



Be s 
+ 9.2818 : 



B 



Bs 



(35) 



again with uncertainties in the coefficients of less than 15 %. The dominant one comes 
from the renormalization scale p\ followed by the CKM factors. For the semileptonic CP 
asymmetries in the Bd system we can also write [46J: 



4 d 



10 4 a 



+ 



Cd + ad 
sin /3 



ad ■ 



cos /3 
+ b d 



- 1 + b d 



cos 2/3 „cos/3 

- 2^r— + 1 



R, 



sin i 



|A„ 



Z' sin 2/3 2 S ^ n ^ 



V ^ 2 



i2* 



cos< 



|A„ 



(36) 



where we have written the (p, rf) dependence in terms of the angle (3 of the unitarity 
triangle and the side R t = ^/(l — p) 2 + rf. 
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The mixing- induced CP asymmetries in Bd — > J/ipKs and B s — > J/ip<f) are very impor- 
tant to constrain (f)^ and S , respectively. For the latter mode an angular analysis is needed 
to separate the different CP components [70]. The mixing- induced CP asymmetries in the 
two modes determine 

sin(<^ + 2/3) and <pf - 2/3 s . (37) 

Here the angle f3 s is defined as positive: 

P s = -arg(-^^J = 0.01818 ±0.00085. (38) 

V * cs * cb J 

(This should be compared with (3 = arg(-V td V tb /V c * d V cb ) « 0.38.) 1 

The measured value |e^ p | implies the following relation among the CKM elements: 



1.25 -10~ 7 = B 



K 



Im [iy cs V* d f A£] Vcc S (^j+lm [(V ts V td f A*] m S (|| 



+ 2Im (V ts V td V cs V; d A%) Vct S 



Here the number on the LHS originates from 



2 2 

m c m t 



M 2 ' M 2 
lvl w w 



(39) 



- 12V2n 2 AM K |e^ p ' 



1.25 -10- 7 = 77^^^, (40) 

GiflmKM 2 , Ke K ' 

The peculiar hierarchy of the CKM elements in Eq. ( 13 9 p enhances the sensitivity to the 
imaginary part of A^. Expanding to lowest non- vanishing order in the Wolfenstein pa- 
rameter A shows 

Im [(V cs V c * d ) 2 A£] = -2A 2 \ 6 f] Re A^ + A 2 ImA£, 

Im [(V ts V* d f Af] = 2A 4 A 10 (1 -p)rj Re A" + A 4 A 10 [(1 -p) 2 - rf ] ImA«, 
21m (y^y^Af) = 2A 2 A 6 r/ReA£ + 2A 2 A 6 (1 - p) Im A£. (41) 

AMk is dominated by physics from low scales. The short-distance contribution is 
dominated by the charm-charm contribution involving the QCD coefficient p cc [61]. There 
is an additional long-distance contribution from box diagrams with two internal up quarks, 
which cannot be calculated reliably. For instance, one could attribute an uncertainty of 
order 100% to the theory prediction of AMk and try to extract a constraint on |A^| 



3 It should be emphasized that the Tevatron experiments, which have presented first determinations 
of <j>^ — 2f3 s from tagged analyses [^nHHO], also use the notation (f> s and /3 S , but with a slightly different 



meaning. Comparing the notation of Ref. 29,30 with our notation one gets 

-2/3 s CDF = tf-2f>.. 



Ref. [29l[30] has neglected 2/3 s in the relation between AT S and <f> s in Eq. (j3]). This is justifed in view of 
the large experimental errors and the smallness of 2j3 s . 
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from AMk- While A c £ is very sensitive to any kind of New Physics which distinguishes 
between the first and second quark generations, we will see in Sect. 12.41 that in MFV 
scenarios all effects on A^ are totally negligible. In an unspecified non-MFV scenario 
both ex and AMk are useless, because A^, A K and A K are uncorrelated with any other 
observable entering the global fit of the unitarity triangle, while in MFV scenarios AMk 
is Standard-Model-like. Therefore we do not consider AMk any further. 

2.4 Three scenarios 

After having introduced our parameterisation of New Physics in terms of the A parame- 
ters in Eqs. (IT7|) and (|39|) . we can now discuss the three different physics scenarios which 
we consider in this article. The common feature of all scenarios is the assumption that 
all relevant effects of New Physics are captured by the A parameters. As long as one 
only considers the quantities entering the global fit of the unitarity triangle in conjunc- 
tion with the observables of B s — B s mixing, this property is fulfilled in many realistic 
extensions of the Standard Model @. However, once a specific model is studied, often other 
quantities (unrelated to the global fit of the unitarity triangle) constrain the parameter 
space; prominent examples are branching ratios of rare decays such as B(B — > X^) and 
B(B S — >■ Such effects cannot be included in a model-independent approach like 

ours. Still, we will see that interesting bounds on the A parameters can be found within 
the broad classes of models defined by our three scenarios. In any specific model covered 
by our scenarios, the constraints on the A parameters can only be stronger, but not weaker 
than those presented in this paper. 

Two scenarios involve the MFV hypothesis. The notion of MFV means that all flavour- 
violation stems from the Yukawa sector. It is usually implied that all flavour- changing 
transitions in the quark sector are solely governed by the CKM matrix, while flavour- 
changing transitions in the lepton sector come with elements of the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix. Strictly speaking, this conclusion is only valid if MFV 
is invoked at or below the GUT scale. If MFV is built into a GUT model at a higher scale, 
it is well possible that imprints of the PMNS matrix can be found in FCNC processes 
of quarks. Indeed, the articles in Ref. [32] discuss supersymmetric GUT models with 
flavour-blind soft SUSY-breaking terms near the Planck scale. The renormalisation group 
evolution involving the large top Yukawa coupling then induces FCNC transitions between 
right-handed bottom and strange quarks at low energies. In our analysis this situation is 
a very special case of the scenario I discussed below. 

2.4.1 Scenario I: Non-MFV 

In this scenario, we do not assume anything about the flavour structure of the New Physics 
interaction. Since here A c £, A K and A K are unrelated to other parameters, we can neither 
derive any constraints on these parameters nor use €k in the global fit. While A^ and 
A s are a-priori independent, the allowed ranges for these parameters are nevertheless 

4 A notable exception are models with large couplings of a light charged-Higgs boson to down-type 
fermion. In such models B(B tv t ), which we assume to be Standard-Model-like, is modified. Another 
exception are models with a non unitary 3x3 CKM matrix, e.g with new fermion generations. 
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correlated through the global fit and the unitarity constraints on the CKM matrix. This 
can be qualitatively understood as follows. Consider a value for |A S | which exhausts the 
range allowed by the hadronic uncertainties in AM S . The good theoretical control over 
the ratio AMd/ AM S then fixes |Vtd| 2 |Arf| quite precisely. The global fit of the unitarity 
triangle further constrains \Vu\ oc y(T— p) 2 + fj 2 , so that a-posteriori the allowed ranges 
for | Ad | and |A a | become correlated. Also the flavour-mixed CP-asymmetry af s measured 
at the TeVatron experiments correlates the parameters A^ and A s . 

2.4.2 Scenario II: MFV with small bottom Yukawa coupling 

We adopt the symmetry-based definition of MFV of Ref. [I] to discuss our two other sce- 
narios. Ignoring the lepton sector here, the starting point is the [£^(3)] 3 flavour symmetry 
of the gauge sector of the Standard Model, which entails the flavour-blindness of this sec- 
tor. The gauge part of the Lagrangian is invariant under independent unitary rotations 
of the left-handed quark doublets Q l L (where i = 1,2,3 labels the generation), and the 
right-handed quark singlets d R and u R in flavour space. In the Standard Model the [{7(3)] 3 
flavour symmetry is broken by the Yukawa interactions. This symmetry breaking permits 
discriminating flavour quantum numbers, quark masses and flavour-changing transitions. 
Within the Standard Model only the top Yukawa coupling y t is large, all other Yukawa 
couplings are small or even tiny. These small parameters pose a challenge to generic exten- 
sions of the Standard Model. This challenge is met by the MFV hypothesis which extends 
the Standard Model assuming that the only sources of [{7(3)] 3 flavour symmetry breaking 
remain the Yukawa couplings. Specifying to the familiar basis of mass eigenstates, we list 
the following consequences of the MFV hypothesis: 

i) Any flavour- changing transition is governed by the same CKM elements as in the 
Standard Model. 

ii) Any chirality flip qn — > qi is proportional to the Yukawa coupling y q (and, by 
hermiticity of the Lagrangian, any q L — > q R flip is proportional to y*). 

iii) Any flavour- changing transition of a right-handed quark involves a factor of the 
corresponding Yukawa coupling. 

iv) FCNC transitions have the same pattern of GIM cancellations as in the Standard 
Model. 

For example, property i) and iii) imply that any bn — > s R transition is of the form 
Vt s V t ly b y*f(\y t \ 2 , \yb\ 2 , \y s \ 2 \), where / is some function of \y t , s ,b\ 2 - Tlie actual power of 
Yukawa couplings in the contribution from a given Feynman diagram is determined by 
the number of chirality flips through property ii). Property iv) ensures that any possi- 
ble contribution proportional to V cs V* b is GIM-suppressed, i.e. proportional to \y c \ 2 , and 
negligible as in the Standard Model. 

However, we deviate in one important aspect from Ref. [I]. We explicitly allow for 
CP- violating phases which do not originate from the Yukawa sector, i.e., we proceed 
as in Ref. [5J. CP violation is an interference phenomenon and involves the differences 
from otherwise unphysical phases. In order to avoid new CP phases one must align the 
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phases of the Yukawa couplings with those of other parameters which are unrelated to 
the Yukawa sector. For instance, in the context of the Minimal Supersymmetric Standard 
Model (MSSM) this tuning of phases affects the \x term, the gaugino mass terms, and 
the trilinear soft SUSY-breaking terms. It is difficult to motivate this alignment from 
symmetries or through a dynamical mechanism. We therefore explicitly permit extra 
CP phases outside the Yukawa sector, i.e. we consider effects from flavour-conserving CP 
phases. Usually such phases are constrained by experimental bounds on electric dipole 
moments, in particular for the MSSM, where stringent bounds on flavour- conserving CP 
phases can only be avoided with quite heavy superparticles. But in our context of generic 
MFV sizeable flavour- conserving CP phases cannot be excluded a priori [25]. 

Refs. [HE] consider two possibilities for the dominant flavour-symmetry breaking mech- 
anism. While the large top Yukawa coupling always breaks the flavour symmetries of the 
gauge sector, one can consider the case that the bottom Yukawa coupling is also large and 
spoils flavour blindness at an equal level (this occurs in the popular MSSM scenarios with 
large tan/3 |71j). Our scenario II corresponds to the case, where only the top Yukawa 
coupling is large. That is, in scenario II we neglect all effects from down-type Yukawa 
couplings. The possible |AU| = 2 operators are discussed in Refs. [HE]. Thanks to MFV 
property iii) four-quark operators with right-handed b or s fields are accompanied by small 
(down-type) Yukawa couplings [U|5] and no other operator than Q in Eq. (jHJ) occurs in 
scenario II. Therefore, the only effect of New Physics is to change the coefficient C in 
Eq. ©. 

An important observation is that C will always be real, even in the presence of flavour- 
conserving CP phases: The MFV hypothesis implies that C is independent of the flavours 
of the external quarks. If we interchange &£ and q^, the corresponding four-quark interac- 
tion will be governed by the same coefficient C and the effective Hamiltonian will contain 
the combination 

{V: q V tb ) 2 Cq Ll ,b L q L ^b L + (V; b V tq ) 2 Cb Ll ,q L b L rq L . 

Now the hermiticity of the Hamiltonian implies C = C*. (An explicit check is provided by 
the MSSM, where flavour-conserving CP-violating parameters enter C only through their 
moduli or real parts). Hence our scenario II corresponds to the case 

A s = A d = A* with </>f = </>£ = </>| A = 0. (42) 

We next discuss an important extension of the MFV analysis of Ref. [I] in the case of 
€k, where we include effects from the charm Yukawa coupling y c . The potential rel- 
evance of these effects becomes clear when one notices that the Inami-Lim functions 
S{m 2 /M^,m 2 t /M^ v ) and Siml/M^y) are proportional to rn 2 /M^. Within the Stan- 
dard Model a substantial contribution to the e^-hyperbola stems from terms which are 
quadratic in y c , and we have to extend our analysis of MFV New Physics to order y 2 . 
Splitting the Wilson coefficients as: 

C = C SM + C NP , C VJ = Cff + Cf, (43) 

MFV constrains the new contributions C NP , Cf^ to obey the following pattern: 



= V u VZV ja V?J{\ yi \*,\vj\*) (44) 



18 



2 Setting the scene 



where / is a real-valued function with f(0,x) = f(x,0) = by the GIM mechanism and, 
of course, C NP (/z) = C^ p (/z). We must now distinguish two cases depending on whether 
the dominant contributions from New Physics affect the diagrams with the light charm or 
up quarks or rather involve particles with a mass similar to or heavier than the top quark. 

In the first case, we have to consider New Physics contributions C^ p and C^ p which 
involve the charm and up quarks on an internal line. Such contributions occur, for exam- 
ple, in box diagrams in which one or both W bosons are replaced by charged Higgs bosons. 
These diagrams with only light internal quark lines lead to negligible effects, if the new 
particle exchanged between the quark lines is a scalar (like a Higgs boson), because scalars 
couple left-handed to right-handed fields and come with the penalty of small Yukawa cou- 
plings. The extra helicity flip (if an internal quark is right-handed) or the GIM mechanism 
(on an internal line with only left-handed quarks) brings in extra Yukawa couplings and 
the contribution to C^ p and C^ p is of order \y c \ 4 or smaller and negligible compared to 
Standard Model contributions, which are proportional to m 2 oc \y c \ 2 - A scaling like in 
the Standard Model, with just two powers of \y c \, could occur in principle if the new ex- 
changed particle is a heavy gauge boson mimicking the Standard Model couplings to the 
left-handed quark doublets. We are not aware of a realistic theory with such particles and 
do not consider this possibility further. 

The second case corresponds to New Physics contributions which involve heavy parti- 
cles and directly add to the coefficient C in Eq. ( ITUl) . This class of contributions includes 
box diagrams with a charged Higgs and one charm and one top quark as internal quark 
lines. These diagrams indeed give a contribution to C^ p proportional to \y c \ 2 - Another 
prominent example for a contribution of this type are the chargino-squark diagrams in 
the MSSM and it is worthwile to discuss this example for illustration, before returning to 
our generic scenario. The chargino-stop box diagram, which contributes to A^, is widely 
discussed in the literature. However, to our knowledge, nobody has studied the corre- 
sponding effect in or A^. The former parameter receives contributions from a box 
diagram with a scharm on one line and a stop on the other. In the limit y c = there 
is an exact GIM cancellation between the contributions from the charm and up squarks, 
which has been invoked to justify the omission of scharm effects. The first non-vanishing 
contribution is proportional to \y c \ 2 (corresponding to cl —> cr and cr — > cl flips). Since 
the same flip is needed on the stop line, there is also a factor of \y t \ 2 involved. Clearly, 
we recognise the pattern of Eq. (jUJ) with the same function / as in A^-. Extending to 
the generic MFV situation, it is easy to relate A^ to A^- for theories in which y t is small 
enough that we can expand A^ and A^ to the lowest order in y t (like in the MSSM 
for moderate values of tan (3 and not-too-large values of the trilinear breaking term A t )^ 
Then 

/(IVilMVif) = fo\Vi\ 2 \Vj\ 2 + 0{\vij\*), (45) 
and up to small corrections one has Cf t p = C^ p ||/ c | 2 /||/t| 2 and C^ p = Cf p \y c \ A /fyt^ ■ Since 
C^ p is real in scenario II, we can certainly neglect it and set A# = 1 in Eq. (1391) . To account 
for the situation that y t is close to one we should vary C^ p around C^ p \y c \ 2 /\y t \ 2 . A realis- 
tic range for C^ p can be obtained from the Standard Model situation. In S(rfi 2 t /M^) the 

5 Thc actual expansion parameter is ytv/M, where M is the mass scale of the new particles in the loop 
and v = 174 GeV is the Higgs vacuum expectation value. 
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relevant quantity is m^ t /M w ~ 2y c>t . However, while S{jn 2 /M^) ~ m 2 c /M^ v , S(fn 2 /M 2 v ) 
differs from m 2 /My V by a bit less than a factor of 2. We take this as a conservative estimate 
and choose 

= \ K ?%p*\c™ with 0.5<A X <2, (46) 

while C^F = 0. As indicated in Eq. (J4"6]) . the Yukawa couplings y c and y t , which are 
expressed in terms of m c and m t , enter the Wilson coefficients at the scale /znp at which 
the heavy particles of the new-physics scenario are integrated out. Since m c /m t is scale- 
independent, we evaluate this ratio at the scale /i = m t in the following. 

Since no new operators occur in scenario II, our New Physics parameters are related 
to the Wilson coefficient C in a simple way: 

r ^NP 

A ij _ _ 1 ( . 7) 



Eqs. P6]) and f )47|) imply the relation 



1 , x " b c\" b t) u \""tl lv± W) 'Itt / Att ^ /, 

1 + Ar- 2 2 „ -,/n,r2\„ l^A'" 1 ) I 4 



mg(m t ) g (m 2 /M^) rju 
m 2 (m t ) S (ml/M 2 ,, rnllM 2 ,) Vct 
1 + 0.017A K (A|-1). (49) 



The numerical value 0.017 is obtained for the central values of Table [71 The smallness 
of this number roots in the enhancement of S '(jn 2 / 'Myy ,m 2 / 'M^) by the large leading 
logarithm log(m^/M^) = —8.3, which stems from box diagrams with internal charm and 
up quarks. This logarithm is absent in the New Physics contribution, which is formally 
of the order of a next-to-leading-order correction. In summary, Eqs. (|42p and (|4"8"|) and 
A c £ = 1 define our scenario II. 

Potential New Physics effects governed by y c have also been studied in [5] in generic 
s — > d transitions. This reference estimates their size as of order 1% of the contribution 
governed by y t and finds them negligible. Our result in Eq. (f49|) . specific to K—K mixing, 
is in agreement with this estimate. In view of the sub-percent experimental error of ex 
and decreasing theoretical uncertainties, the correction A^ is not negligible a priori. 

2.4.3 Scenario III: MFV with a large bottom Yukawa coupling 

In scenario III we consider a large bottom Yukawa coupling. Then h\ AB ^~ 2 in Eq. ((7j) is 
modified to include operators that are not suppressed anymore 

H \AB\=2 = (v ^ )2 + CsQs + + h( , (5Q) 

Here Qs = Ql^rQl^R an d Qs is defined in Eq. f[T3"j) . MFV does not put any constraint 
on Cs and C$ which can be complex (an example of an MSSM scenario with a complex 
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Cg can be found in Ref . [72] , and an up-to-date renormalisation-group analysis of Cs and 
Cs can be found in the appendix of the fifth article in Ref. [16]). In the following, we will 
assume that the matrix elements of Q, Qs, Qs are affected in the same way by [/-spin 
breaking corrections, so that the presence of new operators in ilj/^" 8 ' -2 yields a scenario 
III corresponding to the case A s = A^ with generally non-zero 4>f = <p^. 

Cs and Cs necessarily involve at least two powers of y^ because of MFV property iii). 
In the kaon case, the corresponding contribution to analogous coefficients Cs and Cs in 
jj\i\s\=2 would involve y 2 s instead. Clearly, such contributions to A^- and A^ will have a 
negligible impact on ek, but this need not be the case for A^, in view of the big lever arm 
in Eq. (EH), meaning that the coefficient of ImA^ is larger than that of Re A^ by three 
orders of magnitude. Contributions to A c £ from box diagrams with new heavy particles 
involve four powers of \y c \ in addition to the two powers of y s and are negligible compared 
to the Standard Model contribution oc \y c \ 2 even when multiplied by a factor of 1000. 
New contributions involving internal charm quarks (or some new neutral scalar particle 
coupling a di to an sr) can be proportional to \y c \ 2 as in the Standard Model, but the two 
powers of y s are sufficient to suppress the effect below a level which is relevant for ex- In 
addition, in scenario III, C will be a function of \yb\ 2 , which complicates the relationship 
to the corresponding coefficient in _£/1 AS 1 =2 and even our proof that C is real does not 
hold anymore. In summary, our scenario III comes with A s = A^ and this parameter is 
complex. The parameters A^- are as in scenario II (that is, they are real, fulfill Eq. ( 1481) . 
and A c k = 1), but are now unrelated to A djS . 

We disagree with Ref. [5] on one point here, namely the possibility of (9(50%) effects 
in ex through a sizeable complex contribution to (in our notation) A tt . The claimed effect 
involves two powers of the strange Yukawa coupling y s . The MFV property links this 
contribution to a similar one in B—B mixing, which involves two powers of yi> instead. A 
50% effect in from this source would imply an enhancement of B—B mixing by almost 
two orders of magnitude from this term. We do not see how extra contributions with more 
powers of y^ could possibly reduce this enhancement to a factor below 1. 

One ought to mention that scenarios II and III do not exhaust the possibilities of- 
fered by MFV. For instance, Refs. [THEI] consider MFV-MSSM scenarios with Msusy ^ 
Ma° <; v i where M SU sy and M A o denote the masses of the superparticles and the CP-odd 
Higgs boson, respectively, and v = 174 GeV is the electroweak scale. In the MSSM the 
coefficient Cs is highly suppressed [H] 0, while the operator sl&r sabi, occurs with a size- 
able coefficient, despite of the suppression with the small strange Yukawa coupling. Its 
counterpart in the B d system comes with the even smaller down Yukawa coupling and is 
negligible. In the scenario of Refs. [16j[7l] the connection between B s — B s mixing and 
Bd — B d mixing is lost. But large effects in B s — B s mixing are not allowed, due to the 
experimental bound on B(B S — > /i/z) [16] (for up-to-date results, see [73]). 



6 The vanishing of Cs in the MFV-MSSM scenario with Msusy 3> -M40 > v stems from a softly broken 
Peccei-Quinn symmetry [3], which we have not built into our scenarios II and III. 
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2.4.4 Testing the Standard Model 

There are various ways to test the Standard Model. The simplest one is to determine all 
the relevant parameters from a global fit, and test the fit prediction for a given observable 
compared with the direct measurement. This kind of test, which can be quantified by 
computing the relevant pull value, is independent of any underlying new physics scenario. 

The second kind of test addresses a definite new-physics scenario extending the Stan- 
dard Model, and computes the statistical significance that the parameters take their 
Standard Model value. In our case it corresponds to testing whether one or several A 
parameters are compatible with A = 1. 

In the relevant sections below we will perform both kind of tests, and discuss their 
interpretation. 

3 Inputs 

In this Section, we discuss all relevant experimental and theoretical inputs entering the 
fits. The corresponding values and uncertainties are quoted in Tables [6] and [7J In general, 
if there is only one uncertainty quoted, we understand this error as a statistical one. 
In case of two error contributions, the first one is taken as a statistical error while the 
second, theoretical, error is treated as an allowed range for the observable or the parameter 
under consideration. This kind of uncertainty is treated in the Rfit scheme described in 
Refs. [md. 

3.1 Hadronic parameters and method of averaging 

Several hadronic inputs are required for the fits presented, and we mostly rely on lattice 
QCD simulations to estimate these quantities involving strong interactions at low energies. 
The presence of results from different lattice QCD collaborations with various statistics 
and systematics make it all the more necessary to combine them in a careful way. The 
procedure that we have chosen to determine these lattice averages is as follows: 
We collect the relevant calculations of the quantity that we are interested in and we take 
only unquenched results with 2 or 2+1 dynamical fermions, even those from proceedings 
without a companion article (flagged with a star). In these results, we separate the error 
estimates into a Gaussian part and a flat part that is treated a la Rfit. The Gaussian part 
collects the uncertainties from purely statistical origin, but also the systematics that can 
be controlled and treated in a similar way (e.g., interpolation or fitting in some cases). 
The remaining systematics constitute the Rfit error. If there are several sources of error 
in the Rfit category, we add them linearly 0. The Rfit model is simple but also very 
strict. It amounts to assuming that the theoretical uncertainty is rigorously constrained 
by a mathematical bound that is our only piece of information. If Rfit is taken stricto 
sensu and the individual likelihoods are combined in the usual way (by multiplication), the 
final uncertainty can be underestimated, in particular in the case of marginally compatible 



7 Keeping in mind that in many papers this combination is done in quadrature and the splitting between 
different sources is not published. 
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values. We correct this effect by adopting the following averaging recipe. We first combine 
the Gaussian uncertainties by combining the likelihoods restricted to their Gaussian part. 
Then we assign to this combination the smallest of the individual Rfit uncertainties. The 
underlying idea is twofold: 

• the present state of the art cannot allow us to reach a better theoretical accuracy 
than the best of all estimates. 

• this best estimate should not be penalized by less precise methods (as it would 
happen be the case if one took the dispersion of the individual central values as a 
guess of the combined theoretical uncertainty). 

It should be stressed that the concept of a theoretical uncertainty is ill-defined, and the 
combination of them even more. Thus our approach is only one among the alternatives 
that can be found in the literature [T^JES]- In contrast to some of the latter, ours is 
algorithmic and can be reproduced. Moreover, we differ from the PDG-like method ad- 
vocated in Ref. [76J on two points. We separate systematic and statistical errors, which 
prevents us from assigning a reduced systematics to a combination of several results suffer- 
ing from the same systematic uncertainty. We do not attempt at estimating the (partial) 
correlations between the results from different collaborations, even though we are aware 
of their existence (results from the same gauge configuration, using the same procedure 
to determine the lattice spacing. . . ). Whatever the averaging method chosen, one should 
emphasize that it relies crucially on the quality of the error estimation performed by each 
collaboration. 

The following tables show the inputs used and the average obtained by applying the 
procedure described above for the following hadronic parameters: the decay constant 
fg s for the B s meson (Table [[]), the ratio of decay constants fsj ' fn d (Table [2]), the 
scheme-invariant bag parameter Bb s = 1.523 (m&) for the B s meson discussed af- 
ter Eq. f JT2|) (Table |3]), the ratio of bag parameters BsJBsi (Table @J 0, and the bag 
parameter £>j<-(2GeV) for the neutral kaon (Table [5]). 

We are not aware of lattice estimates for the power-suppressed matrix elements cor- 
responding to <Br ; and Eft.. We will assign them a default value of 1 ± 0.5, taking a flat 
uncertainty for the two bag parameters contributing the most to Ar 9 and a 9 SL (B^ and 
respectively), whereas the remaining ones are either assigned a Gaussian uncertainty 
(Br , Br x , BftJ, or determined from the HQET relations Eq. fTT5]) (£>r 2 , Br 3 ) for which 
we allow a 20% power correction modelized as a flat error. All these bag parameters vary 
independently for Bd and B s mesons, i.e. we do not assume the exact SU(3) symmetry 
for them. 



8 Some of the lattice collaborations [83l[84] provide only the parameter £b = {JbJ fB d )\fBB s /BB d 
together with the ratio of decay constants, without providing the correlation coefficient between the two 
values. In such a case, we have extracted the ratio of bag parameters and its uncertainties errors assuming 
a 100% correlation in the systematic errors between £ B and fB s /fB d - 
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Collaboration 


Nf 


fB a ± O stat ± <?Rfit 


Reference 


CP-PACS01 


2 


242 ± 9±H 


m 


MILC02 


2 


217 ± 6±ij? 


m 


JLQCD03 


2 


215 ± 9t\l 


m 


ETMC09 


2 


243 ± 6 ± 15 


m 


HPQCD03 


2+1 


260 ± 7 ± 39 


m 


FNAL-MILC09 


2+1 


243 ± 6 ± 22 


[82] 


HPQCD09 


2+1 


231 ±5 + 30 


m 


Our average 




231 + 3 + 15 





Table 1: Calculations and average used for the decay constant fs a - Nf stands for the 
number of dynamical flavours used in the simulation. The first uncertainty quotes 
the statistical uncertainty, the second the Rfit error. 



Collaboration 


N f 


fB 3 / fB d ± Vstat ± &Rfit 


Reference 


CP-PACS01 


2 


1.179 + 0.018 + 0.023 


[77] 


MILC02 


2 


1.16 + 0.011^ 


[78] 


JLQCD03 


2 


1.13±0.03^ 2 7 


m 


ETMC09 


2 


1.27 + 0.03 + 0.04 


m 


FNAL-MILC09 


2+1 


1.245 ± 0.028 ± 0.049 


[82] 


HPQCD09 


2+1 


1.226 + 0.020 + 0.033 


[S3] 


RBC/UKQCD10 


2+1 


1.15 + 0.05 + 0.20 


m 


Our average 




1.209 + 0.007 + 0.023 





Table 2: Calculations and average used for the ratio of decay constants fB s /fB d - 
Nf stands for the number of dynamical flavours used in the simulation. The first 
uncertainty quotes the statistical uncertainty, the second the Rfit error. 



Collaboration 


Nf 


Bb s ± V-stat ± ORfit 


Reference 


JLQCD03 


2 


1.299 ± 0.034±g;J|§ 


m 


HPQCD06 


2+1 


1.168 + 0.105 + 0.140 


[85] 


RBC/UKQCD07 


2+1 


1.21 + 0.05 + 0.05 


m 


HPQCD09 


2+1 


1.326 + 0.04 + 0.03 


[83] 


Our average 




1.28 + 0.02 + 0.03 





Table 3: Calculations and average used for the bag parameter &b s - Nf stands for the 
number of dynamical flavours used in the simulation. The first uncertainty quotes 
the statistical uncertainty, the second the Rfit error. 
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Collaboration 


N t 


B B JB Bd ± (Tstat ± CTRfit 


Reference 


JLQCD03 


2 


1.017 ± 0.016t°;°i7 




HPQCD09 


2+1 


1.053 + 0.020 + 0.030 




RBC/UKQCD10 


2+1 


0.96 + 0.02 + 0.03 




Our average 




1.006 + 0.010 + 0.030 





Table 4: Calculations and average used for the bag parameter ratio BbJ^b^ Nf 
stands for the number of dynamical flavours used in the simulation. The first uncer- 
tainty quotes the statistical uncertainty, the second the Rfit error. 



Collaboration 


N f 


Bk(2 GeV) ± a s tat ± o Rfit 


Reference 


JLQCD08 


2 


0.537 ± 0.004 ± 0.072 


m 


HPQCD /UKQCD06 


2+1 


0.618 + 0.018 + 0.179 


m 


RBC/UKQCD07 


2+1 


0.524 + 0.010 + 0.052 


m 


ALVdW09 


2+1 


0.527 + 0.006 + 0.049 


m 


Our average 




0.527 + 0.0031 + 0.049 





Table 5: Calculations and average used for the bag parameter Bk(2 GeV). Nf stands 
for the number of dynamical flavours used in the simulation. The first uncertainty 
quotes the statistical uncertainty, the second the Rfit error. 



3.2 Observables not affected by New Physics in mixing 

In this section, we first discuss the observables allowing us to establish an universal pre- 
ferred region in the p — fj subspace, independent of any New Physics contributions in 
mixing. 

• The CKM matrix element \V U( i\ has been determined from three different methods: 
superallowed nuclear /3-decays, neutron /3-decay and pion /3-decay. Currently, the 
best determination of \V u d\ comes from superallowed /3-decays where the uncertainty 
is dominated by the theoretical error, see e.g. Ref. J9TJ[92]. An analysis by Towner 
and Hardy [93J focusing on an improvement of the isospin-symmetry-breaking terms 
finds a central value for \V ua <\ which is larger though still compatible when compared 
to values quoted in the past, with a slightly reduced uncertainty: \V u d\ = 0.97425 ± 
0.00022. 

• The matrix element \ V US \ can be determined from decays, from hadronic r decays, 
and from semileptonic hyperon decays. We are using the K £3 average quoted by 
Flavianet [91]. The experimental number for f+(q 2 = 0) • \ V US \ obtained by averaging 
results from ISTRA+, KLOE, KTeV, and NA48(/2), as quoted by Flavianet, is 
f + (q 2 = 0) • \V US \ = 0.2163 ± 0.0005 leading to \V US \ = 0.2254 ± 0.0013 [91195]. 

• The matrix element \V c f,\ is obtained from semileptonic decays B — > X c £u, where 
X c is either a D* meson (exclusive method) or a sum over all hadronic final states 



3.2 Observables not affected by New Physics in mixing 



25 



containing charm (inclusive method). For several years the most precise value has 
been provided by the inclusive method where the theoretical uncertainties have been 
pushed below the 2 % level by determining the relevant non-perturbative Heavy 
Quark Expansion (HQE) parameters from moment measurements in B — » X c iv and 
B — > Xs'j decays. The inclusive \V c b\ value used in our analysis is taken from the 
Heavy Flavour Averaging Group (HFAG) [96], |K 6 , incl | = (41.85±0.43±0.59) x 10" 3 , 
where the first error contains the experimental and HQE uncertainties and the second 
reflects the theoretical uncertainty on the total rate prediction for B — > X c iv. 
The theoretical uncertainty on the exclusive \V c b\ determination in the calculation 
of the form factor value at zero recoil F(l) has not been competitive so far. A 
recent calculation provides a significantly smaller error budget: -F(l) = 0.921 ± 
0.013 ± 0.020 [97], although the exclusive \V c b\ determination still gives a larger 
uncertainty. Using the average value for the product F(1)|V^| = (36.04±0.52) x 10 -3 
from HFAG [HS] and applying a 0.7 % QED correction [98], one finds |V^ jCxc i| = 
(38.85 ± 0.56 cxp ± 0.55 thcostat ± 0.84 theosys ) x lO" 3 = (38.85 ± 0.77 ± 0.84) x lO" 3 , 
which has a smaller central value than the inclusive result. 

We average the two \ V c b\ values in such a way that the smallest theoretical uncertainty 
is preserved similarly to our procedure to average lattice inputs and we obtain | V c b\ = 
(40.89±0.38±0.59) x 10~ 3 , keeping in mind, however, that the inclusive and exclusive 
numbers are not in perfect agreement. 

• The two methods to extract \V u b\, the inclusive and the exclusive ones (using the 
theoretically cleanest B — > -niv decays), both suffer from significant theoretical 
uncertainties. The exclusive measurements prefer values around 3.5 x 10~ 3 [96J. 
The numbers quoted are from partial rates measured at large q 2 (> 16 GeV 2 ) or 
at small q 2 (< 16 GeV 2 ), using form factor calculations from lattice QCD [99 |I100] . 
or Light Cone Sum Rules (LCSR) |101j . respectively. The fit input is the average 
of these numbers, following the same procedure as for the lattice QCD parameters: 
iWcil = (3.51 ± 0.10 ± 0.46) x 10~ 3 . 

The average of inclusive results quoted by HFAG using [96] the Shape Function (SF) 
scheme [102] yields (4.32 ±0.16^q 22 ) x 10 -3 , where the first uncertainty contains the 
statistical and experimental systematic uncertainty as well as the modelling errors 
for b — > ulv and b — > civ transitions. Compared to HFAG, we modify the assignment 
of the uncertainties as follows. We add the following uncertainties in quadrature: 
statistical uncertainty, experimental systematics, b — y c and b — )■ u modeling and the 
error from the HQE parameters (6-quark mass nib and /1 2 ). Several theoretical uncer- 
tainties can only be guestimated: the shape function uncertainty, contributions from 
subleading shape functions, weak annihilations and the procedure of scale matching. 
We assign an additional uncertainty on m^, which reflects higher order corrections 
not accounted for in the partial rate predictions for B — > X u tv [103] . We choose 
50 MeV as the additional uncertainty. All these uncertainties of a second type are 
added linearly. As a result, we obtain a significantly larger theoretical uncertainty 
compared to the uncertainty quoted by the HFAG: (4.32l;o; 2 4 ± 0.45) x 10~ 3 . 
The exclusive and the inclusive inputs are then averaged using the same recipe as 
for the lattice QCD parameters and we obtain: \V uh \ = (3.92 ± 0.09 ± 0.45) x 10~ 3 . 
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• A measurement of the branching fraction for B + — > t + v t allows one to constrain 
the product \V u b\ ■ fs where fs is the decay constant of the charged B meson. The 
theoretical prediction for this branching fraction is given by 

B(B ^ T ,) = ^!^l( 1 -J^) 2 |V; t |VIr Bt . (51) 

87T V m B+ J 

We use the experimental value tb+ = 1.639 • 10~ 12 s in our analysis. B{B — > tv) 
combined with the constraint from the oscillation frequency Amj (see Sect. I3.3[) re- 
moves the dependence on the decay constant fs (assuming that the decay constant 
for the charged and neutral B meson is the same, i.e., neglecting isospin-breaking 
effects of order 1%). 

First evidence for the decay B + — > t + v t has been made by the Belle collabora- 
tion 1 104] by reconstructing the decay on the recoil of fully-reconstructed B-meson 
decays. Using the same technique the BABAR collaboration found then a 2.2 a ex- 
cess |105j . By searching for B + — > t + u t on the recoil of semileptonic B-meson 
decays Belle |106] found also evidence for B + — > r + u T while BABAR found a 2.3 a 
excess |T07J. Recently, the Belle analysis first presented in Ref. [106] has been sub- 
mitted for publication [108] reporting a slightly shifted central value and the BABAR 
analysis described in Ref. |105] has been updated to the full data set |105] in which 
evidence has been found for the decay B + — > t + v t . 

The world average for B(B + — > t + u t ) that is used in the analysis calculated from 
the measurements performed by BABAR and Belle [104U10TH109j is B(B + — > t + u t ) = 
(1.68 ± 0.31) x 10" 4 . 

• The input for the CKM angle 7 (= arg[— V ud V* b / iy cd V*^)]) is taken from a combined 
full frequentist analysis of the CKMfitter group using CP-violating asymmetries in 
charged B decays to neutral mesons plus charged mesons. The data are 
taken from HFAG using the three different methods proposed by Gronau, London, 
Wyler (GLW) rjIT], and Atwood, Dunietz, Soni (ADS) [112], and including also the 
Dalitz plot approach developed by Giri, Grossman, Soffer and Zupan (GGSZ), and 
independently by the Belle collaboration [113J. At the 68.3 % confidence level (CL), 
the result of this analysis is (71~^H)° with a second solution at 7 + 71". The constraint 
is shown in Figure [TJ 

The above determination of 7 raises interesting statistical issues. The angle 7 ac- 
tually appears as the complex phase of a suppressed ratio re of decay amplitudes 
(different re's appear for different final states); in other words in the limit tb — > 
there is no constraint left on the CKM phase, and for finite tb the error on 7 is 
roughly inversely proportional to r# itself. It turns out that the current data do not 
exclude tiny values for the re's, one obtains for the DK final state tb = 0.103^qo24 
at 68.3% CL, the D*K final state r B = 0.116^25 at 68 - 3 % CL , and for the final 
state DK* final state r B = 0.111^47 

at 68.3% CL. Because the ratio of amplitudes 
is related to the bare observables non-linearly, its maximum-likelihood estimate is 
biased, and it can be shown that this bias overestimates the value of tb which in 
turn implies an underestimate of the uncertainty on 7. In the statistical language, 
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Full Frequentist treatment on MC basis 



D(*) K(*) GLW + ADS I wa 
DO KO GGSZ □ Combined 
hhkh CKM fit 




Figure 1: Constraint on the angle 7 from a combined analysis of B — > Z)(*).ff(*) 
decays. 



this effect yields a significant undercoverage of the naive 68.3 % CL interval for 7 
computed from the log-likelihood variation. 

A better estimate of the statistical uncertainty on 7 can be obtained by inspecting 
the deviation of the distribution of the log-likelihood among a large number of toy 
experiments from its asymptotic limit. Problems arise because this distribution is 
not only non asymptotic, but also depends on nuisance parameters, that is other 
parameters than 7 that are necessary to compute the toy experiments. In such 
a situation the most conservative approach is called the supremum one, since it 
maximizes the uncertainty over all possible values of the nuisance parameters. To 
date this method which guarantees the coverage properties by construction is the 
default one for the treatment of 7 in CKMfitter, but it must be kept in mind that 
it actually leads to overcoverage in general [114] . 



Another determination of 7 which is unaffected by New Physics in B — B mixing 
is obtained by combining measurements of a and 0: The corresponding quantities 
are changed into a — 0^/2 and (3 + 0^/2 in the presence of a new physics phase 
(f)f which therefore drops out of 7 = 7r — (a — 0^/2) — ((3 + 0^). This procedure 
leads to a significantly more precise determination of 7 than B — >■ D^K^ decays. 
The individual measurements of a — 0^ /2 and (3 + 0^ /2 are described in the next 
section. 
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3.3 Observables in the B^—Bd system affected by New Physics 
in mixing 

The following observables can be affected by |AF| = 2 New Physics contributions in the 
B d ~ B Q d system: 

• The oscillation frequency Am d in the B d sector has been measured with 0(1 %) pre- 
cision mainly due to the 5-factory data [131] . The translation of the measured value 
for Am d into constraints on the CKM parameter combination |V^ d VQ!j| 2 suffers from 
significant uncertainties in the theoretical calculation of the product fs d \fBs d of 
hadronic parameters. These hadronic parameters can be obtained from lattice QCD 
computation, computing fB d \f^B d from f Bs and Bb s , and the flavour-symmetry 
breaking ratios JbJ fB d and BB s /BB d as quoted in Section 13.11 and summarized 
in Table [7J The value and uncertainty for the perturbative QCD correction t)b 
has been originally estimated in Ref. [IS]. With up-to-date values and uncertain- 
ties for a s and the top-quark mass one obtains f] B = 0.8393 ± 0.0034. For the 
top-quark mass we take the average value of the Tevatron Electroweak Working 
Group [115] , m t = (172.4 ± 1.2) GeV, combining published and also preliminary 
results from D0 and CDF. This mass, interpreted as a pole mass, is translated into 
TTLt{jn t ) = (165.017 ± 1.156 ± 0.11) GeV in the MS-scheme at one-loop order. It 
should be noted, however, that the identification of the measured mass value to 
the pole mass is under debate (with potential new systematics coming from this 
identification), see e.g. [116] . 

• In the Standard Model, the predicted decay width difference Ar^ is small: Ar^ = 
(38. ll^) x 10~ 4 ps~ 1 [12J (with our inputs). The average between DELPHI and 
BABAR measurements of the ratio AT d /T d calculated by HFAG [TTTHnglUHT] is 
0.009 ± 0.037. The experimental uncertainty is much larger than the size of the 
Standard Model prediction, AT d /T d = (58^22) x 10~ 4 , and the measured value is in 
good agreement with the Standard Model prediction within experimental uncertain- 
ties so that stringent constraints on New Physics contributions cannot be derived at 
the present stage of precision. Since a huge amount of statistics will be needed to 
measure Ar^ at the level predicted by the Standard Model, this situation will prob- 
ably not change for quite a long time. Even then if a deviation from the Standard 
Model value were observed due to New Physics contributions in mixing it would 
show up beforehand in other observables like ag L or sin 2/3. As a consequence, Ar^ 
has no visible impact in our discussion and it is not used as an input to the fits 
presented here. 

• CP violation in B d mixing (i.e., \q/p\ 7^ 1, with q/p defined in Eq. (JSJ) can be 
measured from the untagged dilepton rate asymmetry 




N e + e + — N1-1 



2(i-k/p|) • 



(52) 



AW + N e -i 
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With a tagged time-dependent decay asymmetry, one measures equivalently 




T(B° q (t)^l+X)-T(B° q (t)^l-X) 
T(B*(t)^l+X) + T(B°(t)^l-X) 

-^ = 2(l-\ q / P \) + 0((\ q /p\-lf) 



(53) 



(54) 



with the time- dependence dropping out. A weighted average of BABAR, Belle and 
CLEO measurements [TT9] results in A§ L = -(47 ± 46) x 10" 4 [T3T] . which is 
a bit less than one standard deviation below the Standard Model prediction of 
aj L = (-7.58^) x 10" 4 POESIES] (with our inputs). 

• Within the Standard Model the measurement of the S = sin 20^ coefficient in 
the time-dependent CP asymmetry A C p{t) = S sin (Am d -t) + C cos (Am d ■ t) in 
decays of neutral B d mesons to final states (cc)K° provides a measurement of the 
parameter sin 2/3, where (3 = arg[— V td V t * b / (V d V* b )], to a very good approximation. 
The current uncertainty of 0.023 on sin 2<\> d is still dominated by statistics [131] . 
The difference between the measured sin 20^ coefficient and sin 2/3 due to penguin 
contributions has been theoretically estimated in Ref. [120] to be below the 10 -3 level, 
while phenomenologically less stringent constraints on this difference are quoted in 

Refs. [rmrzij . 

When interpreting the measured sin 20^ coefficient as sin (2/3 + <ft) the Standard 
Model^FC hypothesis does not rigourously apply. However, the gluonic penguin is 
OZi-suppressed and the Z-penguin is estimated to be small so that New Physics in 
decay is assumed to be negligible with respect to the leading tree amplitude. We 
neglect the effect from possible New Physics in K — K mixing on sin (2/3 + </>^), 
which is justified given the small value of the well-measured CP- violating parameter 

• The measurement of sin 20^* results in two solutions for 2/3 + (in [0,7r]). One 
of these solutions can be excluded by measuring the sign of cos20^ x . For a recent 
review of BABAR and Belle measurements see e.g. Ref. |125j . The current exper- 
imental results from BABAR and Belle using a time- and angular-dependent anal- 
ysis of B° — > J/ipK*° decays, a time-dependent Dalitz-plot analysis in B°/B° — > 
E>(*)o/£,(*)o h o with B o/go ^ D^°/D^°h°, and B°/B° D* + D*~K° S} disfavour 
negative cos 2<p d )K values but HFAG deems it difficult to average the different mea- 
surements or to determine a reliable confidence level as a function of cos20^ x [131] . 
Here, as a simplification, it is only assumed that cos 2<p d hK > 0. 

• The constraint on the CKM angle a = it — (3 — 7 is obtained from time-dependent 
and time-independent measurements in the decays B — > tttt, B — > pp, and B — > pn. 
The time-dependent CP asymmetries measured in B — > tttt provide information 
on the effective parameter sin (2a ef j) (which is a function of a and the penguin- 
to-tree ratio [110] ). It is possible to translate this measurement into a constraint 



30 



3 Inputs 



on a by exploiting isospin symmetry which allows to pin down the penguin-to- 
tree ratio and thus to determine the difference a — a e g from data [126] . Under 
the assumption of exact isospin symmetry the amplitudes A + ~ = A(B° — > 7r + 7r~), 
A 00 = A(B° — > n°n°), and A +0 = A(B + — > satisfy a triangular relationship: 

\/2A +0 — \^2A 00 = A + ~ . A corresponding relationship holds for the CP conjugated 
decays: V2A +0 - V2A 00 = A+~. These isospin triangles can be reconstructed by 
measuring the branching fractions and direct CP asymmetries for the final states 
B° — > 7r + 7r~, B° — > 7i°tt°, and B + — > 7r + 7r°. Since one measures sin (2a e fj) and since 
the triangle has a twofold ambiguity for its apex in the complex plane, there is an 
eightfold ambiguity for a in [0,7r]. The extraction of a from the isospin analysis 
is independent of any possible New Physics contributions in the AI = 1/2 decay 
amplitude except for the singular point a = [127]. If there is no New Physics 
contribution in the AI = 3/2 decay amplitude (as assumed here), the extraction 
provides a = n — 7 — (3 — 4>j /2- As a consequence, a is equivalent to 7 if + (p^ is 
measured e.g. from B — > J/tpK s as already pointed out above. 

Similar in line an isospin analysis can be performed for the B — > pp system. In this 
case the analysis needs to take into account the measured longitudinal polarisation 
of the p mesons in the different final states B° — > p + p~ , B° — > p°p°, and B + — > p + p°. 
Finally, the pir modes provides another crucial input for a, by using a model for the 
Dalitz decay into three pions in addition to the isospin symmetry [128] . The results 
of these analyses that are based on the world averages of the various BABAR and 
Belle measurements for the CP asymmetries and branching fractions determined 
by HFAG [13 lj are displayed in Figure [2j The combined analysis results in a = 
(89.01^)° at 68.3 % CL. ' 

The most stringent constraint on a comes currently from the B° — > pp channel. The 
uncertainty is driven by the rather large branching fraction B(B + — > p + p°). The 
input value for B(B + — > p + p°) has changed recently when the BABAR collaboration 
presented a new analysis on the final data set [129j . The large branching fraction 
value leads to an isospin triangle that just closes. The measurement uncertainty is 
smaller than the expected uncertainty. As a consequence, the current uncertainty 
quoted for a might be on the optimistic side. It should also be stressed that at this 
level of precision so-far neglected uncertainties (electroweak penguins, 7i~— r/ ' mixing, 
p — uj mixing, other isospin violations, finite p-width, etc.) should be considered in 
more detail. 

3.4 Observables in the B s — B s system affected by New Physics 
in mixing 

We now discuss observables which can be possibly affected by |AF| = 2 New Physics 
contributions in the B s -B s system. 

• The experimental input for the mass difference Am s taken from HFAG is dominated 
by the measurement of CDF [2S] • The dependence of the Standard Model prediction 
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Observable 


Value and uncertainties 


Reference 




97425 ± 00022 


[93] 


\V 1 

| v us\ 


2254 ± 0013 




777, ts 


(497.614 ± 0.024) MeV 


[92] 


GrTT 1 
v r 


(1.16637 ±0.00001) x 10~ 5 GeV 2 


[92] 


m w 


(80.398 ± 0.025) GeV 


[92] 




(5.27917 ± 0.00029) GeV 


M\ 
\ i 




(5.27953 ± 0.00033) GeV 


[92] 




(5.3663 ± 0.0006) GeV 


[92] 


\v cb \ 


(40.89 ±0.38 ±0.59) x 10~ 3 


see text 


\v ub \ 


(3.92 ± 0.09 ± 0.45) x 10~ 3 


see text 


B(B -> rv T ) 


(1.68 ±0.31) x 10~ 4 


see text 


/ 


71 +21 ° 

—25 


see text 




(0.507 ± 0.005) ps" 1 


[92] 


air 


(-47 ±46) x 10~ 4 


11311 

i i 




0.673 ± 0.023 


[130] 




positive 


see text 


a 


89.0 ± 4 f 


see text 


Am s 


(17.77 ±0.12) ps- 1 


[28] 




(-85 ± 28) x 10" 4 


see text and [3511132] 




(-17 ±93) x 10~ 4 


[133] 


<$+ vs. Ar s 


see text 


[30ll3T|fT32] 


\£k\ 


(2.229 ±0.010) x IO" 3 


[92] 


Aitik 


(3.483 ±0.006) x lQ- 12 MeV 


[92] 



Table 6: Experimental inputs used in the fits. 
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Theoretical Parameter 


Value and uncertainties 


Reference 


f B 


(231 ± 3 ± 15) MeV 


see Section |3.1| 


Bb (iTlb) 


0.841 ±0.013 ±0.020 


see Section |3.1| 


fBr, 

tl SI t> A - J CL 


1.209 ±0.007 ±0.023 


see Section |3.1| 


B B JB Bi 


1.01 ±0.01 ±0.03 


see Section |3.1| 


Vb 


0.8393 ± 0.0034 


see text 


m t (jn t ) 


(165.017 ± 1.156 ±0.11) GeV 


see text and (115J 


Bk 


(0.724 ±0.004 ±0.067) 


see Section |3.1| 


h 


156.1 MeV 


see text 




0.940 ±0.013 ±0.023 


see text 


Vtt 


0.5765 ±0.0065 


m 


Vet 


0.47 ±0.04 


m 


Vcc 


(1.39 ±0.35) f^GevV 1 

v Trie J 


m 


Xk 


(1.25 ± 0.00 ± 0.75) GeV 


see text 


m c (m c ) 


(1.286 ±0.013 ±0.040) GeV 


see text 


d / OS 


1.01 ±0± 0.03 


11341 

E i 


Bs B. {™>b) 

*-> j J -'S \ u / 


0.91 ±0.03 ±0.12 


|134j 


MS 


(0.222 ± 0.027) GeV 


from a s (M z ) in [92] 


m s (m b ) 


(0.085 ± 0.017) GeV 


m 


m b (m b ) 


(4.248 ± 0.051) GeV 


m 


n?™ 


(4.7±0±0.1)GeV 


m 


Br 


1.0 ±0.5 


m\ 


Br, 


1.0 ±0.5 


m 


B Rl 


1.0 ±0.5 


m 


Br 2 


1.0±0±0.5 


m 


% 


1.0±0±0.5 


m 



Table 7: Theoretical inputs used in the fits. 
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Figure 2: Constraint on the angle a from the isospin analyses of B — > irir and B — >■ pp 
decays and a Dalitz plot analysis of B — > pir decays. In the presence of New Physics 
in B d mixing, the quantity shown is tt — 7 — j3 — 4> d 



Am™ on p — 1] coordinates appears very weak through the relevant CKM matrix 
elements term |V^*V^ 6 | 2 , but the value of Am s gives a direct constraint in the A s 

plane, by computing fs s \J &b s from f Bs and Bb s , given in Table El The hadronic 
matrix element for the B s system can be related to the B d one via the flavour- SU (3) 
breaking correction parameter £ defined through f B BB a = i 2 fB d ^B d Measurements 
of Am s thus reduce the uncertainties on f B Bs d since £ is better known from lattice 
QCD than f B BB d - This relation also generates a strong correlation between the 
New Physics parameters |A^| and |A S | when allowing for New Physics contributions 
to mixing. 

• The Standard Model prediction [12] for a SL is at least one order of magnitude smaller 
than the one for <2g L , see Table [8) Compared to the Standard Model prediction the 
D0 measurement has a quite large uncertainty, a| L = (—17 ± 91I23) x 10~ 4 |133j . 
and hence does not have a strong impact on the New Physics constraints in the B s 
sector. It is nevertheless included in our analysis. 

• D0 |135] using 1 fb" 1 and CDF |132] have measured inclusive dimuon CP asym- 
metries. The D0 result corresponding to the measurement quoted in Ref. |135] 
presented as a^ L + j^a s SL = -0.0028 ± 0.0013 ± 0.0009 is provided in [136] whereas 
CDF [132] quotes the result as 

A SL = f' dZda f + f sZs(lsL = 0.0080 ± 0.0090 ± 0.0068, (55) 

JdZd + J S Z S 

where f d ^ is the fraction of neutral B d ^ mesons in the fragmentation and Z d ^ is 
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given by 

Z » = T^|-TT^ (56) 

with i/g = AT q /2T q and x q = Am q /T q (see also Ref. [31]). Very recently D0 
has presented a new measurement of Asl [35] using 6.1 fb" 1 integrated luminosity 
which shows a 3.2<r deviation from the (almost zero) Standard Model prediction, 
and is the first direct evidence against the Standard Model in B meson observables: 
Asl = —0.00957 ± 0.00251 ± 0.00146. This result supersedes the former result in 
Ref. [136]. The average between the new D0 result and the CDF result reads 

A SL = -0.0085 ± 0.0028, (57) 

which is 2.9 standard deviations away from the Standard Model prediction. 
For the interpretation of the measured observables we use the following values and 
uncertainties for f d ( s \, y q and x q . f d = 0.333 ± 0.030, f s = 0.121 ± 0.015 with 
a correlation coefficient of +0.439 are taken from Ref. [96J, x d is calculated from 
Am d = (0.507 ± 0.005) • lO 12 ^ 1 [92] and r Bd = (1.525 ± 0.009) • 10" 12 s [137], y d is 
set to zero since Ar^ is expected to be very small in the Standard Model and is not 
affected by New Physics in our scenarios, x s is calculated from Am s = (17.77±0.12) • 
lO 12 ^ 1 [92] and y s from t Bs = (1.515 ± 0.034) • 10" 12 s [137], the measurement of 
Ar s in B s — > J/ip (see below) and the flavour-specific B s lifetime t^ s = (1.417 ± 
0.042) ■ 10~ 12 s [137]. This results in Z d = 0.3741 ± 0.0054 and Z s = 1.0044+°;°°| 
hence actually the error on Z djS has negligible impact on our global fits. 

CDF and D0 have presented time-dependent tagged analyses [2"9~|I3"U] of B s — )■ J / 'if) cf> 
decays which provide information on A s through Ar s and 4>f^. These analyses 
supersede the previous untagged studies of Refs. j!38|,ll39j (which use the same data 
sample, but without the tagging information) and we will not use the untagged 
results in our analysis. These time-dependent tagged analyses have raised a lot of 
attention recently, in particular when the UTfit collaboration claimed an evidence of 
New Physics of at least 3a based on a global fit where these measurements played a 
central role |140j . It has been later argued, though, that this conclusion came from 
an overinterpretation of the data [141] [^. In the framework of the HFAG |131j the 
CDF and D0 have determined a combined constraint based on the measurements in 
Refs. [2"9jl3~0] . In 2009 CDF has updated the analysis on a larger dataset using 2.8/6~ 1 
of data [31 J. The new average between D0 [30] and CDF [21] has been presented 
in Summer 2009 [142] . Using this new average the deviation of the measured value 
(f)f^ with respect to the Standard Model value f3 s is essentially unchanged and reads 
2.3 standard deviations. This average is our default input for the corresponding 
observables, supplemented by the constraint on the flavour-specific B s lifetime = 
(1.417 ± 0.042) ■ 10~ 12 [137] which can be viewed as an independent measurement 
of Ar s . New results for B s — > J/ifxf> have been presented in Summer 2010 by CDF 



9 In particular Ref. |140) assumes that the only effect of the SU(3) assumption on the strong phases 
in the D0 analysis amounts to the suppression of a mirror solution, without any impact on the accuracy 
and the location of the main solution. 
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(with 5.2/fe -1 ) [32] and D0 (with 6.1/6 -1 ) [30] collaborations, in closer agreement to 
the Standard Model expectations, but these measurements have not been combined 
together yet. They have not be included in the present analysis, but are briefly 
discussed together with the Standard Model significance tests below. 

3.5 The neutral kaon system 

The master formula for ex has been given in Eq. (139!) . from the relation between ex and 
M( 2 . The translation of ex into a constraint on p and fj suffers from sizeable uncertainties 
in the Wolfenstein parameter A (the determination of which being driven by |V^,| 4 ), Bx 
(see Table [7] and Section I3TT]) . from the long-distance corrections to the relation between 
Mf| and ex encoded in K e , and, though of less importance, from uncertainties in the QCD 
corrections coming from 7] cc [5T], from the charm quark mass m c (m c ) in the MS scheme, 
from m t and the perturbative QCD corrections r]u [EE] and r\ c t [6T] . 

• From the experimental point of view, the number on the LHS of Eq. ( |39|) has shifted 
substantially over time. For instance in 1995 the corresponding number was 1.21 • 
10~ 7 [B2]. More recently, the numerical value for ex has shifted by about 2.3 % 
(a 3.7 a effect) between the 2004 and 2006 edition of the Particle Data Group 
(PDG) from (2.284 ± 0.014) x 10~ 3 [T43] down to (2.232 ± 0.007) x 10~ 3 [92]. This 
shift has been mainly driven by improved measurements of the branching fraction 
B(Kl — > 7r + 7r~) performed by the KTeV, KLOE and NA48 collaborations leading 
to a reduction of 5.5 % of the semileptonic branching fraction values. 

• As discussed in sec. 12. 2\ the relation between ex and K — K mixing is affected by 
several corrections encoded in n e . We have combined them with the experimental 
result for ex as indicated in this section. 

• The kaon decay constant fx is taken from the review on pseudoscalar decay constants 
in Ref. [92J which is calculated from the measured branching fraction B(K + — > 
/i + z/ M (7)) and the measured charged kaon decay time using an external input for 
\V US \. In Ref. [92] \V US \ = 0.2255 ± 0.0019 from decays is used as external input 
which leads to fx = 155.5 ± 0.2 ± 0.8 ± 0.2 MeV where the first error is due to 
the experimental uncertainties, the second due to the uncertainty from \ V US \ and the 
third due to higher order corrections [92]. With our input of |V^ S | = 0.2246 ±0.0012 
this translates into fx = 156.1 ± 0.2 ± 0.6 ± 0.2 MeV. In the fit, we do not consider 
the uncertainties on fx at this point since they are currently negligible with respect 
to the other uncertainties. Since the fx value obtained in this way is anticorrelated 
with our | V us | input a consistent treatment would require including the leptonic kaon 
decay in the fit and constrain fx simultaneously. This would then lead to an even 
smaller uncertainty on fx given the improved uncertainty on \V US \ imposed by the 
global fit. Such a fit is technically possible but has not been performed here since 
the error from fx on ex does not play an important role. 

• Various determinations of the charm quark mass are available. For instance, the 
charm quark mass in the kinetic mass scheme can be obtained from fits to data 
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from lepton energy and hadronic mass moments in B — > X c lv decays combined with 
photon energy moments measured in B — > X s j decays, see e.g. Refs. |131[I144] . 
The most recent value for the kinetic mass quoted by HFAG is m^ m = (1.165 ± 
0.050) GeV |131] . corresponding to a value in the MS scheme: m c (m c ) = (1.265 ± 
0.060±0.050) GeV. The first uncertainty on the charm quark mass is correlated with 
the b-quark mass uncertainty obtained from the same fits quoted in Table |6] with a 
linear correlation coefficient of order 98 %. A second uncertainty of 50 MeV has been 
added following the discussion in Ref. |144j . to take into account the low renormal- 
isation scale and the size of higher-order perturbative corrections when translating 
the mass from one scheme to another. 

As an alternative, the charm quark mass can also be determined from e + e~ an- 
nihilation data into hadrons created from quark-antiquark pairs. The OPE-based 
method consists in writing sum rules for moments of the cross section cr(e + e~ — > cc), 
which are dominated by the perturbative term and the contribution proportional 
to the gluon condensate. An older analysis of Steinhauser and Kuhn based on a 
three-loop calculation finds m c (m c ) = (1.304 ± 0.027) GeV with a small uncer- 
tainty |145j. A similar analysis by Jamin and Hoang |146j obtains a consistent 
result, m c (m c ) = (1.290 ± 0.070) GeV, but quotes a significantly larger uncertainty. 
This can be ascribed to the choice of the OPE scale separating short- and long- 
distance physics and it can be viewed as the impact of (neglected) higher order 
terms in perturbation theory on the determination of m c through such moments. In 
a more recent calculation at four loops for the perturbative contribution, the uncer- 
tainty has been further reduced: m c (m c ) = (1.286 ± 0.013) GeV [HI] B However, 
in this reference, there is only a limited discussion of the freedom in the choice of 
the renormalisation scale for the perturbative series and the gluon condensate is 
varied only in a limited range, even though these two effects were seen as bring- 
ing significant systematics in Ref. |146j . In the absence of further studies on the 
systematics discussed above in the case of the four-loop analysis of Ref. |147] . we 
assign an additional theoretical uncertainty of 0.040 GeV and use as our input value: 
m c (m c ) = (1.286 ± 0.013 ± 0.040) GeV which is consistent with the values from the 
moment fits but has smaller uncertainties. 



4 Quantitative Studies 
4.1 Standard Model fit 

In this section, we present the current status of the Standard Model CKM fit 0. Fig. m 
shows the fit results in the p — fj plane. Table M shows the fit results for various pa- 

10 Using an improved routine for the renormalisation group evolution of the computed value m c (3GeV) = 
(0.986 ± 0.013) GeV the same group hnds a slightly shifted central value m c (m c ) = 1.268 GeV [148] . 

n We stress an important difference in our present definition of the Standard Model CKM fit, with 
respect to more 'traditional' definitions, such as in Ref. [35] : we include here the full input list of Table El 
that is we also take into account B s — > J/4>4> an( A the semileptonic asymetries. The latter observables 
have negligible impact on the determination of the CKM parameters, however, there is some visible effect 
on the fit predictions for the subleading angles and asymetries in the B s system listed in Tables [5] and [5J 



4.1 Standard Model fit 



37 




Figure 3: Constraint on the CKM (p, ff) coordinates from the global Standard Model 
CKM-fit. Regions outside the coloured areas have CL > 95.45 %. For the combined 
fit the yellow area inscribed by the contour line represents points with CL < 95.45 %. 
The shaded area inside this region represents points with CL < 68.3 %. 



rameters. In Tables [8] and [9] we also show the result of the fit for observables that have 
been individually excluded from the fit in order to quantify possible deviations between 
the individual input values and their fit predictions. The good overall agreement of the 
combined Standard Model fit mixes quantities that are in perfect agreement with their 
fit prediction, with others that are individually at odds. Possible deviations between a 
selection of measured observables and their Standard Model predictions are discussed in 
more detail in the following. 

One observes a sizeable discrepancy between the input value of BiB — > tv) (see Ta- 
ble [6]) and its fit prediction (see Table |9]) which is mainly driven by the measured value 
of sin 2(3, and was first discussed in Ref. [10]. Removing either BiB — >■ tv) or sin 2(3 from 
the list of inputs results in a x 2 change that corresponds to 2.9 standard deviations. This 
discrepancy could arise either from a statistical fluctuation in the measured B[B — > tv) 
value, from too small (large) a value of fs d {BB d ), or from New Physics in the B — > tv 
and/or sin 2(3 measurements. There is a specific correlation between sin 2(3 and B[B — > tv) 
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Quantity 


central ± = Icr 


1 r~ij 


1 /^1T 

± UL = oa 


A 
A 

P 
fj 


n si +0.011 

22 5 43 + - 00077 
u.zzo^o_ 00077 

n 144+0-029 

U.l^ -0.018 

0.342±g5ig 


f) 01c; +0.020 
U.010_o.038 

22 54 + - 0015 

U.ZZd^_ 0015 

n 144+0- 054 

U.l^ -0.028 

0-342 lo°; ° 2 3 ° 


n 81 p; +0-029 
U.OIO _o.046 

n 99^,4 +0-0023 
U.ZZd^ -0.0023 

n 144 +0.068 
u - 144 -0.037 

0-342 +£8* 


B K (l) 

Jb s [MeV J {\j 

®b s (!) 
f B Jf Bd (!) 
B B JB Bd (!) 

B*b.W (!) 


0-81 US 

9Q/1 O +9.3 

Z04.y _y g 

1 1 1 cr +0.090 

— U.UoD 

1 91 8 +0-054 
l.Z10_ 046 

1 1Q8+ - 076 
l.lOO _ .094 

-0.6 ±% 


0.81 1^ 

ooir +12 
ZJJ -12 

1.115 +8m 

1 99+0.11 
J-.^^-o.io 

-0.6™ 


0- 811^ 

90c +16 

235 _ 15 
1 1 1 +0.54 

J- • J- 1 n in 
— u. ±u 

1 99+O.I6 

1- 141^ 
-0-6^ 


J [10- 5 ] 


2 98+ - 16 

z - yo -0.18 


9 OS +0.30 


2 Q8 +°- 44 


a [deg] (!) 
/3 [deg] (!) 
7 [deg] (!) 
-2fi s [deg] (! 
-2ft [deg] 


97.6 tH 

28.09 + ?; 4 7 9° 
67.1™ 

) "2.083 Uffi 

9 n8/L+0- 097 
z - uo ^ -0.097 


^— „ _^3 Q 

97.6^o 7 
28.1 tli 
67. lift 

-2.08i8:g 

-2.08 


^ -t-5 Fi 

97.6 ^ 2 5 6 
28.1 + ?i 
67. ll^ 
-2.08^1 

-2.0818:1 



Table 8: Fit results of the Standard Model fit. The notation '(!)' means that the 
fit output represents the indirect constraint, i.e. the corresponding direct input has 
been removed from the analysis. 

in the global fit that is a bit at odds with the direct experimental determination. This is 
best viewed in the (sin 2/3, B[B — > tu)) plane (see Fig. @J, regarding the prediction from 
the global fit without using these measurements. The shape of the correlation can be 
understood more deeply by considering the ratio BiB — > Tu)/Amd, where the decay con- 
stant f Bd cancels, leaving limited theoretical uncertainties restricted to the bag parameter 
B Bd . The formula for the ratio displays explicitly that the correlation between BiB — > tu) 
and the angle <p^ K is controlled by the values of B Bd , and the angles 7 and a = rr — — 7: 

B(B -+ tv) _ 3?r m\ f m 2 T V 1 1 /sin/3\ 2 

Am d ~ 4 m 2 w S{x t ) \ m% + J Tb+ B Bd r] B \V ud \ 2 \ sin 7/ 

The comparison of the indirect prediction of B Bd from the above analytical formula (having 
only B(B —7- tu), Am^, sin 2/3, a, 7 and \V u d\ as inputs, that is an almost completely 
theory-free determination of B Bd ) with the current direct lattice determination B Bd = 
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Figure 4: Constraint in the (sin 2(3, /3(.B — )• tu)) plane. The coloured constraint 
represents the prediction for these quantities from the global fit when these inputs 
are removed while the cross represents the measurements with aid uncertainty. 



1.2691q q9q ^ s gi ven i R Fig- El F° r this test the deviation is 2.9a, dominated first by the 
experimental error on BiB — > tu), a, 7 and second by the theoretical uncertainty on 
B>B d - This tests clearly shows that the semileptonic extraction of \V u b\ has little to do 
with the BiB — > tu) anomaly. Further insight is provided by Fig. where the constraints 



on the decay constant fs d and fB d y ^B d are shown. We compare the fit inputs fs d 



and fB d yBB d taken from LQCD calculations with their predictions from the fit. The 
measured BiB — > tu) value leads to the constraint on fs d represented by the green band. 

The orange band represents the constraint on fB d \j &B d thanks to the Amj measurement. 
The combined prediction for both quantities (red and yellow regions) reveals that the 

predicted value for fs d \J &B d is in very good agreement with the LQCD input. Therefore, 
if the discrepancy is driven by too small a fs d value, the lattice-QCD artefact responsible 
for this underestimation should not affect the more complicated determination of the 

AB 



B 



2 matrix element proportional to fB d yBB d , as already demonstrated in Fig. [5] 
in order to preserve the good agreement between the predicted and calculated values for 



B d - 



Another potential anomaly related to | ex \ has been widely discussed in the literature [8j 
[H], but does not show up with our choice of inputs and statistical treatment. More details 
can be found in the Appendix. Other interesting outcomes of the Standard Model global 
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fit concern the prediction of the recently measured CP-asymmetries by the TeVatron 
experiments, namely in B s — > J/ip(j) and in dimuonic inclusive decays (see Sect. [3]). The 
discrepancy of these measurements with respect to their Standard Model fit prediction, 
together with the B(B — > rv) anomaly, are summarized in Table [HU It is worth noting 
that the Standard Model does not correlate these anomalies between each other, because 
the standard prediction for CP-violation (—20^ and A<$C) in the B s system is essentially 
zero, and hence at leading order has no common parameter with the BiB — > tv) anomaly. 
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Quantity central ± CL = la ± CL = 2a ± CL = 3a 



\v j en 

1 v ud\ {■) 


n 
u 


» 1 -0.00030 


n 
u 


q74 26 + - 00060 

M '^ zu -0.00061 


n 
u 


q7426+ - 00089 

y '^ zu -0.00091 


\Vus\ (!) 





oo ^^+0.00095 
u.uuuyo 





99^4 +0-0019 
u.uuiy 





99^/1+0.0028 

££0^± n nfl9Q 


\V ub \ (!) 





nn^f; +0.00015 

UUOdU _ .00020 





nf) or fi +0.00030 
UUOOU_ .00031 





nn^fi+o- 00046 

uuoou -0.00042 


\Vcd\ 





oo con +0.00077 
zzozy _ Q oo77 





99c; q+0. 0015 
ZZOO_ .0015 





99t:q +0.0023 
z ^< JO -0.0023 


\Vcs\ 





Q7S41 +°-00021 
-0.00018 





Q7341 +0-00039 
-0.00036 





q7 o 41 +0.00057 
y '°^ ± -0.00054 


\V cb \ (!) 





U4£»U8 _ .00528 





+0.0014 
U4di -0.0059 





rwiri +0.0021 
U4Ji -0.0065 


\v td \ 





00861 igs 





00861 + ™ 





nn«fii +0-00044 

UU8Ol_ .00068 


\Vts\ 





04068 





04068 





0407+0-0012 

u ^ u ' -0.0020 


\Vtb\ 





qqqiQc: +0.000057 
yyyioo _o.ooooi8 





qqqi 3 c; +0-000069 
yy»±oo _o.ooo034 





qqqi 3 c; +0-000081 
yyyioo _o.ooO05l 


(pd [deg] 




°-° -3.6 




c o+2.7 
°-°-4.6 




S 8 + 2 - 9 
°-°-5.7 


0s [deg] 


n 422 +0 - 046 

u -^ zz -0.181 


n 422 +°- 098 

U.^ZZ _ 248 


42 + - 16 

u -^ z -0.29 



M [10- 3 ] (!) 
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Table 9: Fit results of the Standard Model fit. The notation '(!)' means that the 
fit output represents the indirect constraint, i.e. the corresponding direct input has 
been removed from the analysis. 
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Figure 5: Comparison between the Bb a parameter used in the global fit and its 
prediction. The coloured constraint represents the prediction of Bs d quantities from 
the global fit while the blue point represents the input value with the corresponding 
1 a uncertainty as used in the global fit. 
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Figure 6: Constraint on fg d and fB d \j&B d - The green band shows the 95.45% 
CL constraint on fs d when B{B — > rv) is included in the fit. The orange band 

represents the 95.45% CL constraint on fs d \J &B d thanks to the Araj measurement. 
The combined 68.3% CL and 95.45% CL constraint is shown in red, respectively, in 

yellow. The black point shows the 1 a uncertainty on fs d and fg d \J B~B d as used in 
the fit. 
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Quantity Deviation 

wrt SM fit wrt Sc. I wrt Sc. II wrt Sc. Ill 
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Table 10: Pull values for selected parameters and observables in the Standard Model 
(SM) and Scenarios I, II, III global fits, in terms of the number of equivalent standard 
deviations between the direct measurement (Tables [6l [7]) and the full indirect fit 
predictions (Tables I8[ 151 [IH H2 | ITU H5| Wf \ [Tgjl . These numbers are computed 
from the x 2 difference with and without the input, interpreted with the appropriate 
number of degrees of freedom. The zero entries are due to the existence of flat 
likelihoods in the Rfit model for theoretical uncertainties |481l74j. The last three 
lines show the pulls for specific correlated combinations of the three most anomalous 
observables. 
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4.2 Scenario I: New Physics in Bd — B^ mixing and B s — B s 
mixing 

In this section, we present the CKM fit for scenario I where New Physics in mixing is 
independently allowed in the B d and B s systems {i.e. A d and A s are independent). These 
fits exclude the constraint from e# since it is not possible to obtain non-trivial constraints 
for the three New Physics parameters in the K sector. The first study of this kind using 
only 5-factory data has been performed in [55J followed by a complete quantitative anal- 
ysis [38] profiting from the large data sets of BABAR and Belle. Analyses taking also into 
account the B s system have been performed by the UTfit collaboration |149|,ll50j . 

In Fig. [7] we show the p — fj plane for this fit, allowing us to constrain the parameters of 
the CKM matrix using observables not affected by New Physics according to our hypoth- 
esis. There are two allowed solutions in p — fj which cannot be distinguished when using 
only \Vud\, \V US \, \V c b\, \V u b\ from semileptonic decays and from B — > tu, 7 and a — (fi d /2 
and (3 + 4> d /2. Once a^ L is added, the second solution at negative p and fj values is clearly 
disfavoured leaving as the only solution the one with positive p and fj values. Figs. [8] 
and [9] show the results of the global fit for scenario I in the complex and A s planes, 
respectively. We emphasize that we assume that and A s are taken as independent 
in this scenario, but that some of the constraints correlate them (such that af s from the 
inclusive dimuon asymmetry, and the ratio Am d /Am s ). Therefore the figures should be 
understood as two-dimensional projections of a single multidimensional fit, and not as 
independent computations. The constraint from Am^ in the ReA^ — ImA^ plane shows 
two allowed ring-like regions. They correspond to the two allowed solutions in the p — fj 
plane when ag L is excluded from the list of inputs. Indeed, in this New Physics scenario, 
Am d is proportional to the product |A d | 2 • |V^V^| 2 , where the second factor is different 
for the two allowed solutions since it is the side of the unitarity triangle connecting (1, 0) 
and {p,fj). The impact of a^ L highlights the power of this measurement to exclude a large 
region of the possible New Physics parameter space even with a measurement precision of 
0(5 ■ 10~ 3 ). In the combined fit the inner ring (which corresponds to the solution for (f)^ 
in the first quadrant near the ImA^ axis) in the complex A^ plane is disfavoured. This 
leaves us with an allowed region for |A^| which is in agreement with the Standard Model 
value A d = 1, albeit with possible deviations up to 40 %. The New Physics phase 4> d , 
mainly driven by the B{B — > tu) vs. sin 2(j)^ K correlation, has the best fit value at —12.9°. 
It can be as large as —21.8° (—27.9°) at the 2 o (3 a) level and shows currently a deviation 
from the Standard Model of 2.7 a. It is interesting to note that the combined individual 
constraint from ag L , a^ L and A^l also favours a negative New Physics phase 4>A d , mainly 
due to the measured negative a^ L value. When B{B — > tu) is excluded from the inputs 
ImAd and hence 0A d is in good agreement with the Standard Model value (see Fig. [TUl) . 
At the same time the allowed range for \A d \ is significantly enlarged since B(B — > tu) 
helps to reduce the uncertainty on Am d : the two rings are enlarged and merge into a 
single one. 

The constraint on |A S | from Am s is more stringent than that for |A^| - thanks to the 
smaller theoretical uncertainty in its prediction compared to Am d - and in good agreement 
with the Standard Model point. It is interesting to note that also for the B s system the 
constraint from B(B — > tu) plays a non-negligible role: when removing this measurement 
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Figure 7: Constraint on the CKM (p, ff) coordinates, using only the inputs which 
are not affected by new physics in mixing: \V u d\, \V US \, \V c b\, \Vub\ from semileptonic 
decays and from B — > tv, 7 (directly and from the combination 7(0) of a — 4>d 1^ 
and f3 + c^/2). Re gions outside the coloured areas have CL > 95.45 %. For the 
combined fit, two solutions are available: the usual solution corresponds to the yellow 
area (points with CL < 95.45 %, the shaded region corresponding to points with 
CL < 68.3 %), and the second solution corresponds to the purple region. 



from the list of inputs the constraint on |A S | becomes weaker since this measurement 
improves the input on the decay constant fs a through the SU(3)-breaking parameter £ 
(compare Fig. [IT] with Fig. [H]). There is evidence for a non-zero New Physics phase <fif at 
the 3.1 a level. This discrepancy is driven by the from D0 and by the <f>f^ analyses 
from Tevatron, but is expected to be somewhat relaxed by the updated measurements of 
[32l[33]. 

We also note that there is an interesting difference in the allowed size of the New Physics 
contribution when comparing the Bd and the B s systems. While in the B s system the size 
of the New Physics contribution is essentially constrained by the Am s measurement alone 
this is not the case in the Bd system. Indeed, the theoretical prediction of Am^ strongly 
depends on the Wolfenstein parameters p and rj whereas this dependence is very weak 
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i — i — i — i — I — i — i — i — r 



excluded area has CL > 0.68 




ReA d 

Figure 8: Constraint on the complex parameter from the fit in scenario I. For 
the individual constraints the coloured areas represent regions with CL < 68.3 %. 
For the combined fit the red area shows the region with CL < 68.3 % while the two 
additional contour lines inscribe the regions with CL < 95.45 %, and CL < 99.73 %, 
respectively. 



for Am s . The constraint on New Physics in i^-mixing relies thus on \V u b\ on one hand 
and 7 on the other hand, the latter being currently dominated by the combination of the 
sin 2(j)^ K and a measurements which is independent of New Physics contributions in B- 
mixing. The theory prediction for the oscillation frequency Am^ depends on the quantity 
I A d | = v /(ReA d ) 2 + (ImArf) 2 . Without a good constraint on |Ad| from other observables 
it can only be predicted with a very large uncertainty as observed in Table [9j The only 
other observables that are sensitive to the modulus of A^ are ag L and Ash but those are 
measured with a precision that is significantly above the Standard Model prediction and 
thus do no constrain very much the range of Am^ (even though they proved powerful in 
eliminating the negative (p,rj) solution) (the same statement holds for Am 8 ). 

Tables [CD and [12] show the fit results for various parameters and observables. We 
also show the result of the fit for quantities that have been individually excluded from 
the fit in order to quantify possible deviations between the individual input values and 
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Figure 9: Constraint on the complex parameter A s from the fit in scenario I. For 
the individual constraints the coloured areas represent regions with CL < 68.3 %. 
For the combined fit the red area shows the region with CL < 68.3 % while the two 
additional contour lines inscribe the regions with CL < 95.45 %, and CL < 99.73 %, 
respectively. 



their fit predictions. The corresponding pull values are listed in Table [101 Among other 
things it is interesting to note that the indirect fit prediction for the dimuonic asymmetry 
Asl = (— 42^g) x 10~ 4 is consistent at 1.2 standard deviations with the D0/CDF average 
(—85 ±28) x 10~ 4 used here, and remains more precise in spite of the uncertainties on the 
theoretical and New Physics parameters. Hence future improvements of this measurement 
are expected to give crucial information on the underlying physics. 

Another important output of our global analysis is the prediction of the difference 
°sl — a SL' that wm be measured by the LHCb experiment in a close future [31]. It 
reads a| L - a^ L = (-39+^) x 1CT 4 (-93 x 1CT 4 < a| L - a^ L < 36 x 1CT 4 at the 3 a 
level), to be compared to the Standard Model result a| L — ag L = (7.93^; 14) x 10 4 
(4.5 x lCT 4 < a| L - ag L < 9.9 x 1(T 4 at 3 a). 

In contrast to the Standard Model fit, our Scenario I relates the and B s anoma- 
lies through the correlated determination of the A parameters. Hence it is particularly 
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interesting to compute the p- values associated with the hypothesis that some specific com- 
bination of the A parameters take their Standard Model value. This is shown in Table [13j 
We have listed several hypotheses because the 'Standard Model' null hypothesis is com- 
posite, i.e. it does not allow to compute the expected distribution of measurements in 
a given set of experiments, because the Standard Model does not predict the value of 
its fundamental parameters. Hence one is a priori free to choose among the numerical 
hypotheses tested in Table [13] the one that models the Standard Model hypothesis. This 
choice introduces some arbitrariness, and thus slightly different answers to the same un- 
derlying question. In the present context, one may view the hypothesis A^ = A s = 1 as 
the most natural choice to represent the Standard Model: this hypothesis is excluded at 
3.6cr by our global analysis. 

However, one should be aware that this hypothesis somewhat dilutes the most anoma- 
lous pieces of information that are related to CP-violation, by including in the test CP- 
conserving components (corresponding to the real part of A^ and A s ). Let us imagine that 
we consider a more general class of models allowing for CP-violation in several different 
processes that are well compatible with the Standard Model. We would have introduced a 
different A parameter for each process. The test corresponding to all A parameters being 
equal to 1 would then receive a small contribution from the three anomalies that we have 
discussed, but this would be hidden by the other measurements in agreement with the 
Standard Model expectations. This illustates why it is sometimes worth testing reduced 
hypotheses, such as Im(Ad) = Im(A s ) = 0, in order to single out specific deviations from 
the Standard Model. The latter CP hypothesis is excluded at 3.8<r by our global analysis. 

We also learn from this table that scenario III to be discussed below, is slightly dis- 
favoured by the data when one considers it as a subcase of scenario I (A d = A s ), in 
agreement with the third paper of Ref. [36J. Finally as already stressed above, the various 
evidences against the Standard Model shown in Table d3] will be relaxed when the new 
Tevatron average of the B s — > J/ip4> tagged analysis is available [32j[33]. However, a very 
rough estimate allows us to predict that at least the Im(Ad) = Im(A s ) = hypothesis {i.e. 
no CP-violating phase in either or B s mixing amplitudes) will remain disfavoured by 
more than three standard deviations. Indeed although the mixing CP-phase is expected 
to come back closer to the Standard Model value [H2JE3], it remains well compatible with 
the indirect constraint from the dimuonic asymmetry, as shown by Fig. [12j In particular, 
the best value for (fif — 2(3 S is only about one standard deviation below the most recent 
CDF and D0 updates for B s — >■ J/ip(j) [32| [33]. not included in the present analysis. Again 
it would be interesting to know the precise form of the combined CDF and D0 likelihoods 
in order to quantify by how much the difference of the indirect global fit prediction of 
with its direct measurement is increased with respect to our estimate of 1.2 a in Table [TUl 
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Table 11: Fit results in the New Physics scenario I. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 
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Table 12: Fit results in the New Physics scenario I. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 
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Hypothesis 


p- value 


Im(A d ) = (ID) 


2.7 


a 


Im(A s ) = (ID) 


3.1 


a 


Ad = 1 (2D) 


2.7 


a 


A s = 1 (2D) 


2.7 


a 


Im(A d ) = Im(A s ) = (2D) 


3.8 


a 


A d = A s (2D) 


2.1 


a 


A d = A s = 1 (4D) 


3.6 


a 



Table 13: p-values for various Standard Model hypotheses in the framework of New 
Physics scenario I, in terms of the number of equivalent standard deviations. These 
numbers are computed from the x 2 difference with and without the hypothesis con- 
straint, interpreted with the appropriate number of degrees of freedom. 
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4 Quantitative Studies 




Figure 10: Fit result for the complex parameter in scenario I when excluding 
B(B — > tv) from the list of inputs. For the combined fit the red area shows the 
region with CL < 68.3 % while the two additional contour lines inscribe the regions 
with CL < 95.45 %, and CL < 99.73 %, respectively. 




Figure 11: Fit result for the complex parameter A s in scenario I when excluding 
B(B — > tv) from the list of inputs. For the combined fit the red area shows the 
region with CL < 68.3 % while the two additional contour lines inscribe the regions 
with CL < 95.45 %, and CL < 99.73 %, respectively. 
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Efl&M —■■ prevASL(WA) 

sulmUiV □ new ASL(WA) < DO 6.1 fb" 1 




<-2P s (deg) 



Figure 12: New Physics scenario I. Indirect constraint on the CP phase 4>f , com- 
pared with the direct measurement by the TeVatron: previous world average |142j 
(in black) and the Summer 2010 CDF and D0 updates [32| l33] (in blue). The dotted 
line represents the full scenario I fit with the previous world average for af s , while the 
green curve is the update after the D0 evidence for a non zero dimuonic asymmetry. 
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4 Quantitative Studies 




Figure 13: Constraint on the real parameter A from the fit in scenario II. The red 
dashed line represents the constraint if only Bj and B s observables are used excluding 
B(B — > tv) and ex- When adding ex the constraint is essentially unchanged (blue 
dotted-dashed line). A significantly stronger constraint is obtained when adding 
B(B — > tv) (dotted green). The full constraint when adding both, B(B — > tv) and 
ex, is shown in green. 



4.3 Scenario II: MFV with small bottom— Yukawa coupling 

In this section, we discuss the MFV scenario which allows to connect the B and kaon 
sectors. Such a kind of numerical analysis has been first presented in [151] . In this 
scenario, there is only one additional real New Physics parameter A, see Eq. (|42|) . As a 
consequence, this scenario has difficulties to describe a situation where the data prefer a 
non-zero New Physics phase in 5-mixing. Indeed the scenario II hypothesis embedded in 
scenario I, that is A d = A s = A with Im(A) = 0, is disfavoured by 3.7 standard deviations. 
The quality of the fit does not change when ex is removed from the list of inputs. 

Fig. [TBI shows the impact of the various constraints on the parameter A. The constraint 
is only slightly changed when adding to the I?-meson observables - where BiB — > tv) 
has been excluded - the ex constraint. On the other hand, when adding B{B — > tv), the 
constraint gets significantly stronger at the la level while at 2a the reduction in the allowed 
parameter space is modest. In Tables [T41 and [TBI we further provide the constraints on 
individual parameters obtained from the combined fit in scenario II as well as predictions 
for the parameters not used as fit inputs. The compatibility of the parameter A with 1 
is good, meaning that this New Physics scenario does not describe the data better than 
the Standard Model, as expected since all discussed anomalies seem to require new CP 
phases. 
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Quantity 


central ±CL = la 


±CL = 2a 


±CL = 3a 


A 
A 

P 

ft 
'/ 


n S1 70 +0.0093 

n 92S43 + - 00077 
u.zzoio_ 000077 

n 145 +°- 0ib 

U-l^O -0.019 

342+^j? 

w,,J ^ t ^' —0.011 


n oi 7+0.019 
U- 81 ' -0.037 

2254 + - 0015 

U.ZZO^ -0.0015 
U- 1 ^ -0.028 
u.<j-*^ —0.024 


n 01 7+0.028 
U-S-L' -0.047 

fl 99 +0-0023 
u.zzu^i _ .0023 

D 14^ +0.066 
U - 140 -0.037 

342 +0 r°i« 

"•"^■^ —0.038 


A 


0-899 IS 


0.90 1^ 


0-901^ 


Bk 0) 

f Bs [MeV] (!) 
Bb. (!) 

fB s /fB d (!) 

B B JB Bd (!) 
Bs, Bs {m b ) (!) 


0.881^| 
254.5 r^ 4 9 

u - ' 00 -0.075 

1 91 7+0.053 
L -^ L ' -0.035 

1 1Qfi +0-076 
i ' iou -0.095 

-0.9 1£2 


0.8818:1* 

254 1* 7 
0-761°:?° 

1 91 7+O.II3 
-0.076 

1.141^ 
-0.93: 6 7 


0.881^ 
2541^ 

76 + 1 - 12 

u - ' u -0.19 

1 99+O.I6 
ilZZ -0.18 

1.141^ 
-0.9 It? 


J [io- 5 ] 


2 qq+ - 15 

-0.12 


2 qq+ - 29 


2.99 1^ 


oi [deg] (!) 
/3 [deg] (!) 
7 [deg] (!) 


93.3 t 5 d 
28.053; 2 7 ? 
67.1 1 2 . J 


93.3l|? 
28.1 + _H 

67.i^;l 


93.3™ 
28.1 l?i 

67 1 +6 -° 
°' - 1 -10.0 



0rf [deg] 




c 7+2-5 


c 7+2.8 
— °- ' -5.2 


-2A [deg] (!) 


Z.UOO _ .095 


-2.083H 


-2.08 +_° f 7 


-2/3 s [deg] 


9 ns^ +0-070 

Z.UOO _ .097 


-2.08 18:^ 


-2.081°:^ 


<Ps [deg] 


n 4m +°- 042 
_ 124 


0.401 l™?? 


0-4018;^ 



Table 14: Fit results in the New Physics scenario II. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 
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Quantity 



central ±CL = la ±CL = 2a 



±CL = 3a 



\V ud 

\v us 
\v ub 

\Vcd\ 
\Vcs\ 

\v*\(i) 

\Vtd\ 

\Vts\ 

\v tb \ 



U - y ' 4 ^ b -0.00030 

fl OOK/lc; +0.00095 
U.z.z.ou -0.00095 

U.UUOO I _ .00014 

n 22529 +a00077 
u.zzozy _ 00077 

97341 + - 00018 
U - 3,M1 -0.00018 



0.0493 



+0.0028 
-0.0061 



0.00863 lo™ 



0.0407; 



, +0.00038 
' -0.00097 



0.999131 



+0.000040 
-0.000016 



0.974261™ 
0.2254 iKg 

nns57 +0 - 00030 

U.UUOO I _o. 00029 

fl OOKQ +0.0015 
u ' zzdo -0.0015 

+0.00036 
-0.00036 



0.97341 



0493 +0 0048 
u.u^yo_ .oo94 

0.00 8 63 + a00032 
0.04078 



0.999131 



0.00055 

0.00076 
0.00171 

+0.000071 
-0.000032 



0.97426 

fl 99^4 +0-0028 
U.ZZO^_ .0029 

nnS57+ 00046 
U.UUOO / _o.00044 

fl OO +0.0023 
u.z,z,oo _o.0023 

97341 + - 00054 

U.»(O £ ±l_ .00054 



0.0493 



+0.0066 
-0.0116 



o.oo863 ±8:888^ 



0.040£ 



-0.0011 
-0.0021 



0.999131 



+0.000088 
-0.000048 



6K [10- 3 ] 

Am d [ps~ 
Am s [ps~ 

a sl [io- 4 ; 
a S l [io- 4 ; 

a? 
a: 



SL 


— a 


SL 


d 

SL 


[10 


-4] 


d 

SL 


[10 


- 4 : 


s 

SL 


[10 


- 4 ; 


s 

SL 


[io- 


- 4 ] 



(!) 
] (!) 
] (!) 
(!) 

[IO" 4 ] 
0) 

0) 



Ar d [ ps 
Ar s [ps- 
Ar s [ ps - 



-ii 



(0 



i 07+0.54 
-0.55 

n prpr 4 +0.073 
U.OO^_o.o47 



16.2 



-1.7 
-1.4 



a n4 +i.oi 

4 - U4 -0.53 

-4 06 +0 96 

^• uu -0.55 

o 74 +0.99 
°- ' 4 -1.97 



-8.36 
-8.37 
0.373 



+1.93 
-0.94 

+1.90 
-0.95 

+0.045 
-0.083 



fl Q7Q +0.045 
u -°' -0.083 

00426 + 00056 

U.UU £ ±ZU_ . 00154 

O.llO^ 
0.0946 loS 



1 07+0.95 

J-- 8 ' -0.67 



0.554 



1-0.114 
-0.095 



16.2 i 3 ; 6 7 



-4.0 



+1.9 
-1.1 



_4 1+1-9 

-1.1 



1.9 
4.0 

+3.9 
-1.8 

+3.8 
-1.8 

+0.078 
-0.168 



8.7 
-8.4 
-8.4 
0.373 



n Q7Q +0.078 
U.O I O _q.168 

0043 +0 0011 

U.UUOO _ 0021 



0.110 



h0.089 
-0.037 



0-095 l™ 47 



1 07+1.21 

u - dd -0.14 



16.2 



H8.5 
-3.6 



-4-ol?:| 
-4.H?:J 



+4.7 
-2.7 

+4.7 
-2.7 

+0.11 
-0.21 

+0.11 
-0.21 



-8.4 
-8.4 
0.37 
0.37 



nnzn+ - 0016 

U.UUOO _ 0024 

o.nolo ;^ 
o-095±8:~? 



B(B -+ tv) [10- 4 ] (!) 
B(B ->• rz/) [IO" 4 ] 



653 + - 277 

U.UOO_ 00 40 

92 +0 - 12 
u - y ^ -0.10 



fi^+ - 404 

U.UOO_ .077 

92 +0 - 23 

u - yz ' -0.26 



0.65 



+0.52 
-0.11 



0.92 t° il 



Table 15: Fit results in the New Physics scenario II. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 
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Hypothesis 


p- value 


Im(A) = (ID) 


3.5 a 


A = 1 (2D) 


3.3 a 



Table 16: p-values for various Standard Model hypotheses in the framework of New 
Physics scenario III, in terms of the number of equivalent standard deviations. These 
numbers are computed from the x 2 difference with and without the hypothesis con- 
straint, interpreted with the appropriate number of degrees of freedom. 

4.4 Scenario III: Generic MFV 

In Fig. [14] we present the result from the combined fit to Bd and B s observables in the 
complex plane A = A s = A^ for the MFV hypothesis allowing for a large bottom Yukawa 
coupling yt> (scenario III). In scenario III, A can have a sizeable complex component 
proportional to y\. One expects that the constraint A s = Ad reduces the size of the 
allowed New Physics contributions significantly with respect to the general case studied 
in scenario I. However, our fit result in Fig. still allows for a New Physics contribution 
of order —20 % to +40 %. The New Physics phase in this scenario shows evidence for a 
deviation from the Standard Model since both af s and the 4>f measurement on the one 
hand and the discrepancy between sin 2<fi d )K and B(B — > tv) on the other hand, prefer a 
negative New Physics phase in Bd, s mixing. 

Compared to scenario II this model is in much better agreement with the data regard- 
less of whether one takes into account B(B — > tv) as an input or not. While the Bd sector 
constraints do not allow too large a New Physics phase, the B s sector prefers a large New 
Physics phase, though with large errors. As a consequence and as already discussed above, 
the quality of the fit in scenario III is manifestly worse (2.1 a) than in scenario I. 

In this scenario, the Bd- and 5 s -systems decouple from the kaon system and hence 
the constraint on ReA and ImA cannot be improved by adding ex as this introduces 
the additional NP parameter A K in the theoretical prediction for ex- However, when 
including ex in the fit one is able to constrain A K from the fitted values of the CKM 
parameters. The parameter is in agreement with the Standard Model expectation but 
could differ by more than 40 % from A K = 1. Similarly to the previous scenarios, one- 
dimensional results from the fit in scenario III are shown in Tables [17] and [18] Again 
one sees that in particular the semileptonic asymmetries ag L and Ash are more precisely 
predicted than the measurements, so that improvements of the data will be extremely 
instructive. The asymmetry difference is now predicted to be a§ L — a§ L = (18.6 ^i5° 2 ) X 10~ 4 
(— 1 x 10~ 4 < cig L — cig L < 25 x 10 -4 at 3 a), hence a moderate positive value is preferred 
with respect to Scenario I that predicts a larger negative value. 

Finally, two tests of the Standard Model are interesting to study within scenario III, 
and are shown in Table [16] As in scenario I, they show evidence for New Physics, and 
that scenario III describes the data significantly better than either the Standard Model or 
scenario II, assuming as before that the improved data on B s — > J/ip4> do not change the 
overall picture dramatically. 



4 Quantitative Studies 




Figure 14: Constraint on the complex parameter A = A^ = A s from the fit in 
scenario III. For the individual constraints the coloured areas represent regions with 
CL < 68.3 %. In addition, a CL < 95.45 % contour is shown for the individual 
constraints obtained from Am^ and Am s , from A$l-, a sL> ano - a SL> ana - f rom 4>d K 
and (f)f . For the combined fit the red area shows the region with CL < 68.3 % 
while the two additional contour lines inscribe the regions with CL < 95.45 % and 
CL < 99.73 %, respectively. 
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Quantity 




central ±CL = la 


±CL = 2cr 


±CL = 3a 


A 
A 

P 
V 




7928 +0 - 0288 
u. ( yzo _o.0092 

U.ZZOOU_ oo77 

n i 77+0.032 
"J- 1 ' ' -0.012 

0.439 lo :^ 


0-7931^ 
0.2254 lo ^ 5 

n i 77+O.O59 
u - ± ' ' -0.026 

0-439 S 


n 7QQ +0-050 

u. i yo _ .027 
fl 99 +0-0023 

0.1771°:™ 
0.439 Iffi 


Re(A) 
Im(A) 
|A| 

i A r t i 

A [deg] 

K 




0.876 IS 
-0 226 +0 ' 057 

U.ZZU „ .048 

Qn7+ - 128 
u.yui _ .070 

-1 A A 4-9 Q 

-14.4 I^i 

i.oi+H 

—0.38 


0-87618:$ 
-0 23+ 011 

U-^ -0.12 

n q 1+0.35 

u - yi -0.11 

-14.4 Zi\l 

—0.43 


0-88 ±g;g 
-0.23^:11 
91 +0 - 45 

U - yi -0.14 

* ^ 4-1 1 d 

-14.4 1^ 4 

—0.47 


fa. [MeV] (!) 

Bb. (!) 
jBs/ JB d y-j 
BbJB Ba (!) 

Bs, Bs {m b ) (!) 




274 1 2 ^ 

1 44+0-33 
J-- 44 -0.61 

1 230 +0 - 066 

i i +0.090 
l.lOO _ .138 

1 14+0-98 


274 1 4 ? 

1 44+0-99 
-■-•^^ -1.37 

1 23 +ai3 

1 - Z " J -0.10 
i - iO -0.24 
J-. 1 -4.0 


274l 8 7 3 7 

i-4li:? 

1 23 +0 - 18 

- L - Z " J -0.21 

1 16+ - 28 

J-- 10 -0.40 

1 1 +4 - 6 
J-- 1 -5.4 


J [io- 5 ] 




o fi i +0.16 
-0.14 


o fil +0.27 
°- ui -0.37 


3-61 lo :^ 


a [deg] (!) 
(3 [deg] (!) 
7 [deg] (!) 




88.0 HI 

28.01±?;|? 

68.111J 


88.0±g : § 
28.01^ 
68.1™ 


88.0 1^ 4 

28.0 1 2 : ° 
68.1 1JI 


A + 2/3 [deg] (!) 
<Pd [deg] 

A - 2/3, [deg] (!) 
A - 2& [deg] 
0s [deg] 


32.7+I;t 

-20.9 ±2:1 

-16.61L 7 
-17.1+1? 
-13.911? 


qq+16 
°°-29 

-20.9 tH 
-16.6l 8 i 

-17.1 Its 
-13.9 1^ 


qq+27 
oo _ 45 

-2111? 

-16.6™ 

-17 1 +1L8 

1 f .1 _g 4 

-13.9™ 



Table 17: Fit results in the New Physics scenario III. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 
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Quantity central ±CL = la ±CL = 2a ±CL = 3a 



\v ud \(\) 





074 Qn +o-ooo30 
y / <±ou _ 00030 





074^0 +0-00060 
y i _ 00060 





07430 +0-00° 89 

y IIOU _0.00090 


\Vus\(\) 





oocrq/f +0.00095 
ZZJJ ^-0. 00095 





99c;q +0.0019 
z,zoo _o.0019 





99CO +0.0028 
ZZOO _o.0029 


\v ub \(\) 





nn^77+°-ooo 45 

UU<J I I -0.00068 





nnc,77+o-ooo 86 

UUOf t _o.00124 





oo58ig:881? 


\Vcd\ 





99 co/l +0.00077 
ZZOZ £ ±_ .00077 





99^9+0.0015 
z,z,jz, _ .ooi5 





99^9 +0.0023 
ZZOZ _ .0023 


\Vcs\ 





07347 +0-00018 
y '°^' -0.00018 





07047 +0.00036 
y -0.00037 





07^/1 7 +0.00053 
y ' °^ ' -0.00056 







0090 +0.0036 
UOZO _0.0035 





HQ9Q +0.0096 
UOZO_ 0076 





n o 9 +0.017 
uoz -0.011 


\Vtd\ 





nns/i 7 +0.00032 

UU84 ( „o.00032 





00847 +0-00046 
UU84 ( _ 0006I 





OO8/I7+0-00 058 
UU84 f _ .00076 


\Vts\ 





03QR1 +0 00137 
uoyu±_ . 00037 





OQOfil +0.00184 
UOyUl _o.00074 





oqqfi +0.0022 


\v*\ 





000170+0-000016 

yyyi <y_o.oooo57 





000 1 70 +0.000031 

yyyi 1 y _ 000077 





000170+0.000046 

yyyi 1 y _o. 000094 



Am d [ps- 1 ] (!) 


S62 +0 - 081 

U.OUZ _ .068 


o.se 18:11 


0.56 +°io 8 


Am s [ps" 1 ] (!) 


14.91?;? 


14.9 1 4 ;* 


14.9 1™ 


Asl [10- 4 ] (!) 


-22.6 l 3 7 i 


-22.6 1 12 9 9 


90+19 

zo -14 


A S l [10- 4 ] 


-23.4 ±1* 


-23.4 1 1 ^ 


-231H 




18.6l 2 ° 2 


18.6 1 4 ^ 


18.6l? 9 3 4 


4l [10- 4 ] (!) 


-32.4 t£? 


-32.4 1 15 7 9 


09 +22 
JZ -12 


<4 [10- 4 ] 


09 pr+10.2 


09 r +15.5 

OZj.O y y 


00 +22 
00 _ 12 


«Il [10" 4 ](!) 


-13.9+2;! 


— 13.9 I13 2 


_14+ 12 
-■-+-19 


«sl [ID" 4 ] 


-13.9 H 4 


— 13.9 1 13 2 


_14+ 12 
ift -19 


Ar d [ps- 1 ] (!) 


0.00432 lo o o°?8 5 o 3 


0043 +0 0012 

u - uu ^° -0.0023 


43 +0- 0015 

U.UU^9_ .0027 


Ar s [ps- 1 ] (!) 


164+0- 021 

u - lu + -0.074 


lfi4+0- 031 

U.lU £ ±_o,097 


164 + - 040 

U.1U £ ±_ .109 


Ar s [ps- 1 ] 


1 00 +0- 023 

U.1UU_ 019 


100+0- 072 

U.1UU_ 34 


mo+o- 090 

U.lUU_oo46 



B(B^ru) [10- 4 ] (!) 1.36l °i° 1-3618; 
B(B ->• rz/) [10- 4 ] 1.381°;^ 1.381°; 

Table 18: Fit results in the New Physics scenario III. The notation '(!)' means that 
the fit output represents the indirect constraint, i.e. the corresponding direct input 
has been removed from the analysis. 



1.36 



1-0.31 
-0.77 



1 00 +0.31 
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5 Conclusions 

In this paper, we have studied three different scenarios of New Physics in |AF| = 2 
transitions. The complex parameters quantifying the New Physics contributions to B q —B q 
mixing are A q = \A q \ ■ e l ^ d = Mf 2 /M 12 ,q - In K — K mixing three parameters A K , A K , 
and A^ are needed. 

We have first recalled the result of the Standard Model fit using the current available 
data sets. In the B system an interesting effect is observed: the inclusion of BiB — > 
tv) obtained from a combination of BABAR and Belle's measurements deviates from its 
prediction in the Standard Model fit by 2.9 cr which either points to a large statistical 
fluctuation in the experimental numbers, to a problem in the calculation of both the decay 
constant fs d and the bag parameter Bs d on the Lattice, or to New Physics contributions 
in sin 2/3 and/or in B — > tv. If there were New Physics contributions to B& mixing this 
discrepancy would point to a negative no n- vanishing New Physics phase 0^. A second 
hint for a deviation is observed in the B s — B s mixing phase 4>f — 2(3 S measured in B s — > 
J/ip4> though the significance is here only around 2.3 a. The largest discrepancy actually 
comes from the recent improved measurement of the dimuonic asymmetry by the D0 
collaboration, which is at odds at the 3.2 a level with respect to the indirect fit prediction 
(2.9 a when the average with the CDF measurement of the same quantity is made). 
Furthermore, the correction factor n t [9] in the theoretical prediction of ex decreases the 
quality of the Standard Model fit. However, with our inputs (especially the range for Bk 
in Tab. [7]) and with the conservative Rfit error treatment of theoretical uncertainties used 
in our fit we do not observe a significant deviation between the measured ex value and its 
prediction from a Standard Model fit excluding the ex measurement. 

In our New Physics scenario I, we have considered uncorrelated New Physics contribu- 
tions to Bd, B s and K mixing. That is, the complex parameters A^ and A s are allowed 
to vary independently and the kaon sector is omitted, since the New Physics parameters 
A K , A K , and A^ are unrelated to all other observables entering the fit. The experimen- 
tal data are well described in this scenario which can accomodate negative New Physics 
phases preferred by a) the discrepancy between B(B — > tv) and sin 2/3 both measured at 
the B-factories BABAR and Belle, b) the 2/3 s measurements in B s — > J/ip<f) a t the Tevatron, 
and c) the dimuon asymmetry af s measured by D0. The size of the New Physics contribu- 
tion both in B s -mixing and B^-mixing can be as large as 40 % and the hypothesis of zero 
new CP phases, <fid = <Ps = 0, is disfavoured by as much as 3.8 standard deviations in this 
scenario (see Tab. [TBI . The large parameter region emphasises that, despite of the success 
of the B-factories and the Tevatron, there is still considerable room for New Physics in 
Bd as well as in B s mixing. 

In addition, we have considered two Minimal Flavour Violation scenarios. The first 
MFV scenario, scenario II, corresponds to small bottom Yukawa couplings, leading to 
Ad = A s = A K = A with all New Physics phases identical to zero, = </>^ = 4> K A = 0. 
This scenario A^ = A s = A K = A with vanishing New Physics phases has been widely 
studied in the literature. In this scenario, the New Physics parameter A^ is equal to one 
to a very good approximation and A K only deviates by a few percent from one where 
the deviation can be estimated in terms of A K . Since in this scenario no New Physics 
phases are allowed, the deviations seen in the Standard Model fit are still present. As 
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a consequence, this MFV scenario is currently disfavoured by 3.7 standard deviations, 
but not totally excluded yet. The New Physics parameter A can deviate from one by 
about +40 % at 95 % C.L.. The constraint gets only slighly relaxed when removing either 
B(B — > tu) or €k from the inputs to the fit. 

Our scenario III is a generic MFV scenario with large bottom Yukawa coupling and 
arbitrary new flavour-blind CP phases. In this scenario, the kaon sector is unrelated to 
the .B-sector. As in scenario II, one has A^ = A s = A, however, this time with a complex 
parameter A. The New Physics parameters in ex are unrelated to A and can be 
removed from the input list. However, when including ex one is able to constrain the New 
Physics parameter A K which is found to be consistent with the Standard Model value of 
1, but can deviate from unity by about ±40 %. The fit describes the data significantly 
better than the Standard Model fit and than the fit in scenario II since the data prefer a 
negative New Physics phase in Bj and in B s mixings. The hypothesis of a zero new CP 
phase, Im A = 0, is disfavoured by 3.5 standard deviations (see Tab. [T6|) . As in the other 
New Physics scenarios, the allowed size of the New Physics contribution in B^ and in B s 
mixing can be as large as +40 %. 

The several scenarios discussed here show that we have indeed sensitivity to New 
Physics in the |AF| = 2 sector. While the overall picture of current data reveals strong 
hints for New Physics, the current experimental uncertainties prevent us from excluding 
the Standard Model, as highlighted by the p-values of each hypothesis. It has to be seen 
how this picture evolves with the improvement of both experimental and theoretical re- 
sults. When we completed this study, new results from the Tevatron experiments were 
given for the measurement of (f)f^, in better agreement with the Standard Model, which 
will be included in our analysis once the CDF and D0 collaborations have agreed on a 
combination of their results. Importantly, more precise measurement of the CP asymme- 
tries in flavour-specific decays, a SL , from either the Tevatron or the LHC experiments, may 
become crucial in the future: for the time being the theory-and-data driven fit predictions 
for cig L are more precise than the direct measurements, as can be verified by comparing e.g. 
Tabs. |6]and[T2j Hence future more precise data on ag L could help to discriminate between 
the Standard Model and different scenarios of New Physics. Meanwhile, our predictions 
of ag L in Tabs. [T21 dH and [T7] are an important side result of our analyses: they permit a 
fast extraction of a| L from future measurements of <2f s (or a SL — ag L considered by LHCb), 
which is more accurate than what is obtained using the experimental value of af s listed in 
Tab. El 

As an illustration of this statement, we show the indirect fit prediction for the difference 
function of (fif — 2(3 S , in the Standard Model and the New Physics Scenarios 
I and III in Fig. [T5J The prediction in Scenario I and III differs by about two standard 
deviations, hence a precise direct measurement of either observable could not only exclude 
the Standard Model, but also select one of the New Physics scenarios. Thanks to the 
specific two-dimensional correlation it could also invalidate all scenarios, which would 
then imply that there are other sources of New Physics than just contributions to the 
mixing amplitudes. The LHCb experiment is expected to measure both CP asymmetries, 
with an accuracy of about cx(0f — 2(3 S ) = 4° and cr(a SL — ag L ) = 20 x 10~ 4 for 1 fb _1 of 
integrated luminosity [34] . 
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Figure 15: Indirect fit prediction for the difference a s sl — a d sl as a function of 4>^ — 2/3 s , 
in the Standard Model and the New Physics Scenarios I and III. The allowed regions 
correspond to CL < 95.45 %. The prediction from Scenario II is not shown since it 
is very close to the SM. Note that for this plot, the direct measurement of cj)^ — 2f3 s 
in B s — > J/tp(p is removed from the inputs. 



On the other hand, significant improvements on the measurement of BiB — > tv) can 
only be expected from a Super- B-like electron-positron machine since the BABAR and Belle 
results do already rely on most of the available statistics collected at the B-factories PEP- 
II and KEKb. Since the theoretical translation of the ratio BiB — > Tu)/AMd into a 
constraint on CKM and New Physics parameters relies on the calculation of the decay 
constant fs d and the bag parameter Bb 4 , future progress in Lattice QCD calculations is 
also very important. The same remark applies to the hadronic matrix elements entering 
M( 2 and r* 2 . We hope that the current exciting experimental situation will stimulate 
novel activities in lattice gauge theory in this direction. 



64 



A Relationship between ex and K — K mixing 



ACKNOWLEDGMENTS 

We thank D. Becirevic, D. Guadagnoli, L. Lellouch, P. Marquard, M. Neubert and N. Uralt- 
sev for stimulating discussions and comments. We are grateful to A. Kagan for bringing 
a mistake in Eq. (j4"5]l in an earlier version of this paper to our attention. We would also 
like to thank the members of the CKMfitter group for their interest and support for the 
various aspects of this article. 

This work was partially supported by the ANR contract ANR-06-JCJC-0056 and the 
EU Contract No. MRTN-CT-2006-035482, FLAVIAnet. The work of UN is supported 
by the DFG through project C6 of the collective research centre SFB-TR9 and by BMBF 
grant 05H09VKF. 



A Relationship between €k and K — K mixing 

A.l Corrections to the usual ek formula 

In Sect. 12. 2\ we have discussed the connection between ex and K—K mixing. The resulting 
equation, Eq. ( 12 6 p . has been obtained thanks to several approximations (concerning the 
phase 4> e , neglect of £, computation of M12 from lowest-dimension operators of the effective 
Hamiltonian) . The corrections to these approximations are encoded in the deviation of 
the factor n e from 1. A series of papers [9~ll59~ ll6"6T, l6"T] has assessed more precisely the value 
of this factor in order to account for the terms neglected by the previous approximations. 

f ^ and (j) e ^ 45° imply K e ^ 1, Ref. [9J finds K e = 0.92 ± 0.02. In Ref. [76] a 
lattice QCD calculation of I111/I2 is combined with the experimental value of e' K /eK to 
compute £ and finds agreement with Ref. [9] while quoting an even smaller uncertainty: 
n e = 0.92 ± 0.01. Finally, correcting for ImM^ 7^ ImMjj by including higher-order 
terms of the operator product expansion leads to n t = 0.94 ± 0.02. Actually, in the 
correction factor K e , the three approximations have to be treated in different ways, since 
they mix uncertainties from experimental and theoretical origins. The correction from <j) £ 
is of experimental nature and can be treated easily. 

The correction involving £ could be computed by combining the experimental value of 
Re Aq and the theoretical computation of Im A using the effective Hamiltonian if l AS 1 =1 . 
However, the latter is dominated by the matrix element of the QCD penguin in the 1 = 
channel ((7T7t)o|Q6|-^)> which is poorly known. Here we follow the method of Refs. [9l67p76j 
which uses e' K to correlate ((tt^oIQqIK) with A2. The latter amplitude involves the matrix 
element of the electroweak penguin Q$ in the 1 = 2 channel, ((im)2\Qs\K) , which has been 
computed using lattice simulations and sum rules. Indeed, one finds [9] 

— = — -(1 - ft), with u = — — , Q = -T—7-- (59) 

ex V2\e K \ ReA ulmA 

uj = 0.0450 is obtained from experiment, whereas i? describes the weight of the (imaginary 
part of the) 1 = 2 contribution with respect to the 1 = one. 

In practice, a numerical estimate of the contributions from other (subleading) opera- 
tors than Qq and Qs in ijl AS l =1 has been obtained [65lll52Lfl53j under the assumptions 
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that Im Aq and Im A 2 can be computed accurately combining the effective Hamiltonian 
approach and experimental values for Re A Q and Re A 2 . 

e J£ = No + Nl . R( , + N2 . R&: R& = R 6 [(e'/e) cxp , R s ] (60) 

where N , Ni, N 2 are numbers, coming mainly from A t , the experimental values for the real 
parts of A and A 2 , and R$ and i?§ are rescaled bag parameters. Following the review of 
ref. [154] , the authors of ref. |152j proposed the following conservative range Rs = 1.0 ±0.2 
which we will follow. Using e' K /eK = (1.65 ±0.26) • 10~ 3 , one can use eq. (l6"Ui) to determine 
Q, which corresponds to the ratio between the 1 = 2 and 1 = contributions to e' K /ex- 
In principle, this would require to split iVo into two pieces coming respectively from 1 = 
and 1 = 2, and to assess the size of these contributions following Ref. [65]. A quicker way 
to obtain a similar estimate consists in computing |153j : 

where Q\ and Q 2 correspond to the (extreme) hypothesis that N Q is saturated either by / = 
or / = 2 contributions. Eq. fl59|) can be then used to compute £ in either hypotheses, and 
we will take the spread of the obtained values as a (conservative) systematic uncertainty 
in the determination of £. A more detailed analysis of the contributions to iVo would allow 
us to decrease this systematic uncertainty. 

The last correction comes from higher-dimension contributions to ImM 12 . As shown in 
ref. [66], there are two different corrections at d = 8, corresponding to the AS = 2 d = 8 
operators and the double insertion of AS* = 1 operators connected by a u, c loop. The 
first is expected to be very suppressed compared to the d = 6 contributions, whereas the 
second one is essentially dominated by long-distance pion exchanges estimated in Chiral 
Perturbation Theory for weak processes, leading to: 



e K = sin </> e e^ e 



ImM 



(6) 
12 



AM 



p = 0.6 ±0.3. (62) 



For our purposes, it will be simpler to combine these corrections with the experimental 
input for e K : 

M _ ( 1 J K 1 ^ _ e K,exp 



K VV2sin0 £ e K u{l-Q)J n t 

where denote the approximate value of e# in eq. ( )26|) with k 6 = 1. Depending on the 
choice of i? (i.e. the respective part of I = and 1 = 2 contributions in the formula for 

e'/e), we get 

«W = 0.943 ± 0.003(e'/e) ± O.O12(0 e ) ± 0.004(i? 8 ) ±0.015(p) (64) 

nf ] = 0.940 ±O.OO3(e7e)±O.O12(0 £ )±O.OO4(i2 8 )±O.O18(p) (65) 

Combining the first two errors in quadrature for the Gaussian part and the last two errors 
linearly and taking the spread of the values into the Rfit part, we obtain the estimate 

n e = 0.940 ±0.013 ±0.023 (66) 

in good agreement with n e = 0.94 ± 0.02 in ref. [66J. This determination relies on the 
assumption that e' K is not affected by New Physics. 
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Figure 16: Constraint on \ex\ predicted from the global fit, compared to the experi- 
mental value. The green curve is obtained with our Rfit treatment of the uncertainties 
for the parameters entering Eq. (|26f) . The dashed curve is obtained by treating all 
the errors as Gaussian. 



A. 2 Error budget for ex 

Recently, it has been claimed that there is a discrepancy between the ex constraint and 
its prediction 019]. With the input values used in our fit and with the Rfit treatment 
of theoretical uncertainties we do not observe any sizeable discrepancy when comparing 
the prediction for ex (Table E]) with its measured value (Table E]). This can also be seen 
in Fig. [3] where the combined fit prefers a region in p and rj that is close to the edge of, 
though still inside the 95 % CL region of the ex constraint. 

As illustrated in Fig. [16] one can obtain a minor discrepancy at 1.2 a if one treats as 
Gaussian all the parameters involved in Eq. ( l26l) (i.e., \V c b\, Bx, K e , the QCD correction 
terms r] cc ^ tt and the masses m c ,i), but such a treatment of the systematics remains ques- 
tionable. Another way of seeing the absence of discrepancy is to compare the prediction 
for hi t and estimates of these quantity. One can see clearly from Fig. [ITJthat the global fit 
can cope easily with a value of n e down to 0.9, and that the prediction from the global fit 
agrees well with the recent estimates of this quantity. 

In order to make the comparison of our \ex\ analysis with the one of Ref. [76] easier 
we treat all the errors as Gaussian and calculate the error budget for the fit prediction of 
\ex\- With our inputs we find 

10 3 |e^| = 1.893 ±0.020| A ±0.020| A ±0.063| (M) 
±0.180| Bjr ± 0.019|to P ± 0.084| 

charm ±0.054L 

= 1-89^:1 ' ' (67) 




which is 1.2 a away from the experimental measurement, while with the inputs of the last 
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Figure 17: Constraint on K e predicted from the global fit using Eq. ()26|) . compared 
to three different theoretical inputs n e = 0.94 ± 0.02 with either a Rfit or a Gaus- 
sian treatment of the error, or K e = 0.94 ± 0.013 ± 0.023 corresponding to our own 
treatment of the uncertainties involved in its evaluation. 



reference in [76J we find 

10 3 |e x | = 1.769 ±0.019U±0.019| A ±0.061| (M) 
±0.062| Ba . ± 0.018| top ± 0.067| 

charm 

±0.033L 

= l-77l8:il (68) 

which is 2.4 a away from the experimental measurement, in agreement with [75]. Hence 
we see that the difference between our analysis and Ref. [75] mainly comes from our input 
for 13 k that has a larger theoretical error, and from our central value for |e#| that is 
larger because of the different analysis of the CKM parameters. We thus conclude that 
the potential anomaly in |e#| cannot yet be precisely quantified independently of the 
theoretical inputs and therefore deserves further investigations. 
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